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ANNALEN DER PHYSIK 


7, FOLGE * BAND5, HEFT 7-8 « 1960 


Statistische Formulierung von chemischen Reaktionen 
zwischen Fehlstellenteilchen in kubischen Kristallgittern 


Von Alfred Scholz 


Inhaltsübersicht 


Die Gleichungen von Stärstellenreaktionen Gin. (19) und 
(26)) in kubischen Kristallen werden mit Hilfe von Wahrscheinlichkeitsan- 
sätzen (s. Gln. (2) und (9)) abgeleitet. Hierbei fungiert als entscheidende Größe 
die Aufenthaltswahrscheinlichkeit eines Störatoms auf einem Gitterplatz; 
für diese wird eine statistische Deutung (s. Gl. (1)) gegeben. Bei den chemischen 
Reaktionen zwischen Störstellen wird die Korrelation nächst benachbarter 
Fremdionen berücksichtigt. 


Konkret durchgerechnet werden: 
1. die Bildung von Lückenpaaren aus isolierten Lücken in Kristallen vom 
NaCl-Typ, 
2. die chemische Reaktion Ags + Se = [AögSe] in Silberhalogeniden. 
Die statistische Ableitung ergibt für die kinetischen Gleichungen eine ge- 
wisse Modifikation gegenüber dem sonst üblichen phänomenologischen An- 


satz für die kinetischen Gleichungen. Dies betrifft die kinetischen Konstan- 
ten, welche noch eine Zeitabhängigkeit (s. Gln. (19) und (26)) enthalten. 


I. Einleitung 


Bei der atomistischen Berechnung der Beweglichkeit und Diffusions- 
konstanten von Störatomen in Kristallgittern beschränkt man sich im all- 
gemeinen auf den eindimensionalen Fall!). Die entscheidende Größe ist dann 
die Anzahl N(x) der Fehlstellen in einer Gitterebene mit der Ortskoordinate 
x?). Für Diffusions- und Leitfähigkeitsprobleme reicht dieses Verfahren völlig 
aus. Dagegen bei der Berechnung chemischer Reaktionen zwischen verschie- 
denen Störstellenteilchen muß man dreidimensional rechnen. Dann kann man 
nicht mehr die Anzahl N (r) der Fehlstellen in einem Volumenelement AV als 
wesentliche Größe verwenden, weil N (r) eine kleine Zahl wird. Die Volumen- 
elemente müssen nämlich im dreidimensionalen Fall hinreichend klein gewählt 
werden infolge starker räumlicher Änderung der Zahlen N (tr). 

Statt dessen fungiert nun als entscheidende Größe die Wahrscheinlich- 
keit f(r,t) dafür, daß auf einem Gitterplatz (G.P.) mit dem Ortsvektor r 
zur Zeit t ein Störteilchen sich befindet. 


1) Siehe hierzu beispielsweise: A. B. Liddiard, Philos. Mag. 46, 1218 (1955). 
2) N (x) muß eine große Zahl sein, damit die Streuungen klein werden. 
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Die Wahrscheinlichkeit f (r, ¢) kann statistisch so sehn werden: 
Es sei eine gedachte (virtuelle) statistische Gesamtheit, bestehend aus N 
(mit N > 1) identisch aufgebauten Kristallen, gegeben. Dann gilt offenbar: 


N (tt) 
fet) = im, () 
Hierbei bedeutet 


N (t, t) = Anzahl der Kristallgitter innerhalb der gegebenen statistischen 
Gesamtheit, bei denen alle zur Zeit t auf dem G.P. r ein Störatom sich befindet, 

Ähnlich können auch die im Abschnitt II eingeführten Bedingungswahr. 
scheinlichkeiten P (t,, t,, £) statistisch erklärt werden. Mit Hilfe dieser stati- 
stischen Interpretation können die Ausgangsgleichungen (2), (9), (14) abge- 
leitet werden. 

Die Bewegung der Störteilchen im Kristall beruht auf Sprüngen der 
Fehlteilchen von einem Gitterplatz zu einem der bunnehkarten Gitterplätze 
(bzw. Zwischengitterplätze). 


Konkret berechnet werden in dieser Arbeit: 


1. Die chemische Reaktion der Bildung von Läckenpaaren ı aus isolierten 
Lücken. 


2. Die Erzeugung von Assoziaten [AgöSg] in Silberhalogeniden aus 
Silberionen auf Zwischengitterplatz (Ago) und auf G.P. eingebauten $;- 
Fremdionen. 


aay = Durch die Rechnungen werden abgeleitet: 

1. Die kinetischen Reaktionsgleichungen, wobei gewisse Änderungen auf- 
treten, die in den phänomenologischen Ansätzen dafür nicht berücksichtigt 
sind. 

Rich 2. Die kinetischen Konstanten, welche noch eine Zeitabhängigkeit zeigen; 


wt oo 3. Die Diffusionskonstanten von isolierten Lücken, Liickenpaare und 


Ago-lIonen. 


II. Wahrscheinlichkeitsansitze für die chemischen Reaktionen 


Sämtliche Rechnungen werden für einfache kubische Kristallgitter (vom 
NaCl-Typ) durchgeführt. Hierzu werden kartesische Koordinatensysteme 
{x, y,z} mit den Koordinatenachsen parallel zu den Kristallachsen benutzt. 
Bei solch einem kartesischen Koordinatensystem seien 1,, i,, i, die Einheits- 
vektoren, welche in x-, y-, z-Richtung zeigen. Der Abstand # zwischen 
nächstbenachbarten G.P. wird im verwendeten Koordinatensystem stets als 
Längeneinheit angesehen: 


Längeneinheit 


& 


d — Gitterkonstante. 


a) Die Wahrscheinlichkeitsansätze für chemische Reaktionen 
werden am Beispiel der Bildung von Lückenpaaren behandelt. Vorgegeben 
sei ein einfach kubisches Gitter mit A+-Ionen auf den +-Gitterplätzen und 


B--Ionen auf den —-Gitterplätzen. Z.B. kann At -- Nat und B-+- CT] 4, 
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einander entsprechen. In diesem Kristallgitter AG Bg seien nun Ay- und Bp- 
Lücken (Gitterleerstellen) vorhanden. Zur Zeit t = 0 seien alle Ay- und Bp- 
Lücken isoliert. Zur Zeit ¢ = 0 mögen die isolierten A, auf wenigen dis- 
kreten +-G.P. mit den Ortsvektoren r, und die isolierten BJ auf —-G.P. 
mit den Ortsvektoren r;, sich befinden. Im Laufe der Zeit entstehen Lük- 
kenpaare. Es gibt nun 6 verschiedene Orientierungen im Raume für die er- 
zeugten Lückenpaare. 


Für die 
maßgebend: 
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weiteren Abhandlungen sind folgende Wahrscheinlichkeiten 


fa (ut) = Wahrscheinlichkeit dafür, daß zur Zeit ¢ eine isolierte Ay-Lücke 
auf dem -+-G.P. mit Ortsvektor r = {z, y, z} sich befindet; 


fa (t', t) = Wahrscheinlichkeit dafür, daß zur Zeit ¢ eine isolierte By-Liicke 
auf dem —-G.P. mit Ortsvektor tr’ = {2’, y’, 2’} sich aufhält; 
y; (t,t) = Wahrscheinlichkeit dafür, daß zur Zeit t ein Lückenpaar der Orien- 


tierung j existiert mit der assoziierten Ap, auf dem -+--G.P. x sowie 


mit der assoziierten Bp auf dem —-G.P. (t + ;,); 
i, = Einheitsvektor der VORNE A im Raum; 
i; hat die Werte +7,, +i,. +7, 


b) In nullter Näherung wird folgender MER EEE 
benutzt: 
(t, = ha (t, t) - Ip (t (2) 
3): 
Py (tp (0), t + (3) 
u 


Es bestehen die Beziehungen 


= und fp (t 


(r, (0), r, £) Fin): 
Ty 


P , (t, (0), r, £) ist die Bedingungswahrscheinlichkeit dafür, daß eine Ay-Lücke 
zur Zeit t auf dem +-G.P. r sich aufhält, wobei vorausgesetzt wird, daß diese 
Ap-Liicke zur Zeit 0 auf dem +-G.P. r, sich er: 


Das Ziel der Rechnung ist die Bestimmung von a 7 (tr, t). Dazu benötigt 


man zunächst: 


+ At): fp lt + ty, 6+, At) — fz (t + (4) 
Nun gelten aber die Gleichungen: 
f, (t,t-+ At) = = fa (to, t) - (to v, At) + fy (t, Py 


Pp (ta (t), t At) + fp (t+ 
tlt + th) 
»P,(t + + Ab). 3 (5) 


®) In diesen Gleichungen wurde vorausgesetzt, daß im Laufe der Zeit keine neuen 
A, bzw. B, entstehen oder RR (z. B.an Ve rsetzungen). 
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Für kleine Zeitintervalle At haben die Bedingungswahrscheinlichkeiten P 
folgende Gestalt): 


P 4 (to (t), t, At) = & (tp, t)- At de, — + 0(AP) 
P, (t (t), t, At) = 1 — At + 0 (AP) 
Pp (ta (t), + tg, At) = B (to, ty) At + 0 (At?) 


Pp (t + ig t + ig, At) = 1 — Byes At + 0 (AP). 


(6) 


_ Hierbei ist: a (t,, t) - At = Wahrscheinlichkeit dafür, daß die A-; durch 
einen einzigen Sprung während der Zeit At gelangt vom G.P. rt, nach dem be- 
nachbarten G.P. r, wobei ja |r, — r| = ) 2 gelten muß. 

Die Ausdriicke (5) werden in die Formel (4) eingesetzt. Dazu benutzt man 
fiir die Bedingungswahrscheinlichkeiten noch (6). AnschlieBend geht man zur 
Grenze At > 0 über und erhält so: 


é 
—w, (rv, t) = f, (v, t) - Bi, c+%)- — 
a Wx ( ) ha ) fp ( 8) B (te 1 x) (t+ ty 
fp (t + fa (Up, (Tp, - 1,2 — Ve (t, £) - [Xges + Bges] 

nn Im ersten Term von (7) kommen auch Werte r= r+ i; mit i;+%, 
vor. Solche Werte stellen Lückenpaare y, (rt, t) dar. Analoges gilt auch für 
den zweiten Term. Im Ansatz (2) wurde die Korrelation zwischen nächstbe- 


nachbarten Lücken A-, und Bp (Lückenpaar!) nicht berücksichtigt. 


Weiterhin kommt im Ansatz (2) nicht klar der Unterschied zwischen iso- 
lierten und assoziierten Liicken zum Ausdruck. 


Man kann versuchen, die Korrelation zwischen nächstbenachbarten A, und B, 
(Liickenpaar!) durch einen modifizierten Produktansatz zu berücksichtigen: 


ty 


P% (tv, t,t) ist hierbei eine Bedingungswahrscheinlichkeit, welche angibt: Die isolierte 
A, soll vom G.P.r, (zur Zeit 0) zum G.P. r (zur Zeit t) gelangen, vorausgesetzt, daß einer 
der nächst benachbarten —-G.P. (t + i,) von einer B, besetzt ist. 
P* hat eine entsprechende Bedeutung. 
Leider führt der Ansatz (8) zu einem falschen Resultat für + y, (t,t). Der Grund 
hierfür ist: Im Term f, (tr, t) - ¥ p* (tp, t + i,t) muß sich die A, zeitlich zuerst auf dem 
r’ 
u 
+-G.P. r isoliert befinden. Die B, muß danach zu einem späteren Zeitpunkt noch hinzu- 
kommen auf den —-G.P. (rt + i,). Diese zeitliche Nacheinanderfolge ist aber in der Formel 
(8) nicht berücksichtigt. Es gelingt also nicht, yp; (t, £) durch einen modifizierten Produkt- 
ansatz (wie z. B. in (8)) darzustellen. 
c) Bei der richtigen ersten Näherung wird die Korrelation zwischen 
nächstbenachbarten Lücken Ag und Bp berücksichtigt, indem man solch ein 
Lückenpaar [A- BJ] als ein neues Störteilchen y, (t, ¢) darstellt. 


*) Oc, — ns2ist Kroneckersymbol: also für (tg — t)? = 2ist 6 = 1 und für (1 — tr)? +2 
ist 6 = 0. 0 (At?) ist Ordnungssymbol: Es gibt an, daß die weiteren Glieder von Größen- 
ordnung At? und höher sind. 
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Der Wahrscheinlichkeitsansatz der ersten Näherung lautet nun?): 


[y; + At) — (r, t)] fs (v + t;, t) 5 fa (To, t)- (to; t, At) + 


(t+ ty — ix) 


SS fg +i, A+ S y(t, +i, 
ip) ig Fi; 
t+i,A)+ w(t +i,—i,,6)- B, (e+ AN — 
(9) 
—y, (r, t)- sx P44 (t,t, At)—y; (t,t): SS +i,, Yo, At) — 
— (t, 8): P, (t,t + At) — 


j 


7 (r, £) P, (r + 0 + lt) 


P¥ (tp, t, At), Py (tp, t, At), Pi (tr, At) stellen stets Bedingungswahr- 
scheinlichkeiten dar, für die Wanderung einer Ag vom +-G.P. r, (zur Zeit t) 
nach dem +-G.P. r während der Zeit At. Und zwar muß bei P*, Ap auf 
+-G.P. r, isoliert, auf +-G.P. r assoziiert sein; bei P, muß A, assoziiert 
sein auf beiden G.P. r, und r. Bei P muß Ap auf +-G.P. r, assoziiert und 
auf G.P. x isoliert sein. 

Der Wahrscheinlichkeitsansatz (9) kann mit Hilfe der Vorstellung über die 
statistische Gesamtheit (Gl. (1)) het werden. 


Aus der Gl. (9) erhält man nun für z a Vi (rt, £) den Ausdruck®) 


0 
a = (t +t, > =, (To, ale +t; Or, — + 


4 r, t) > to, t)- —t)- r, t)- B* 

+ ip + y,(t + F t) - + iz — 

: ce Ber (10) 

y; ( ) (| (t+ tj > nor 


-y,(r,t)- =. 0,2 — Y; (r, - ß*- + i; 
Fi; 


5) In Gl. (9) soll-+--- andeuten, daß noch weitere Glieder vorkommen müssen. 
Solche weggelassenen Glieder sind aber alle von Größenordnung (At)? und höher, haben 


deshalb keinen Einfluß auf Sy; (r,t). Z.B. wurde nicht berücksichtigt der Fall: Zur 


Zeit t ist weder A, auf +-G.P. noch B, auf —-G.P. (rt + i,). Beide Lücken gelangen 
während der Zeit At auf diese G.P. r und (t + ij). 

6) Die Bedeutung der Sprungfrequenzen a, ß, «*, ß* usw. ist aus der Tab. 2 zu er- 
chen. Bei der Berechnung von (10) werden "für die Bedingungswahrscheinlichkeiten. 


Pi (t, to, At) usw. ähnliche Beziehungen wie in (6) benutzt. Rear Baar 
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In den Formeln (9) u. (10) wurde nicht berücksichtigt: Ein Lückenpaar [4,B,] kann 
noch weiterhin entstehen, indem zu einer isolierten A, (bzw. B,) eine ursprünglich asso. 
ziierte B, (bzw. A,) hinspringt. Ferner kann ein Lückenpaar entstehen, indem zu einer 
isolierten A, eine B, hinspringt, welche ursprünglich zu einem höheren Liickenaggregat?) 
gehörte. Die umgekehrten Vorgänge für die Vernichtung eines Lückenpaares sind auch 
möglich. Solche Fälle haben nur sehr geringe Wahrscheinlichkeit, weshalb sie vernach- 
lässigt werden. 

In der ersten Näherung werden isolierte Aj-und By-Lücken, deren Abstände 
> 1 sind, als unabhängig voneinander angesehen. Die Korrelation zwischen 
ihnen wird vernachlässigt. In der Gl. (10) sind dies die Ag und By mit dem 


Abstand )/3 bzw. //5. Will man auch zwischen ihnen eine Korrelation be- 
rücksichtigen, so muß man beide Lücken als einen angeregten Zustand 
eines Lückenpaares auffassen. Dies ist Aufgabe der höheren Näherungen. 


d) In der zweiten Näherung wird so die Korrelation zwischen A; 
und Bs mit der Entfernung V3 berücksichtigt. Die beiden Lücken Ag und 
B- mit dem Abstand J 3 stellen ein Lückenpaar im 1. angeregten Zustand dar. 


kommt dies zum Ausdruck, indem f, (rt. ¢) - fp, (8). 
wobei |t/, = Y 3 sein muß, ersetzt wird durch x (t, t). Dabei soll x (r, £) be- 
bedeuten die Wahrscheinlichkeit für die Existenz eines Liickenpaares im 1. 
angeregten Zustand mit der Ag auf dem +-G.P. r zur Zeit t. 


Jedoch werden die höheren Näherungen nicht detailliert behandelt. 


Die Wahrscheinlichkeitsansätze wurden am Beispiel der Lückenpaar- 
bildung dargestellt. Aber auch für andere chemische Reaktionen von Fehl. 
stellen in Kristallen gelten analoge Relationen. 


4 III. Kinetische Gleichungen für die Entstehung von Lückenpaaren 


In einem einfachen kubischen Gitter A&Bz gibt es sechs verschiedene räum- 
liche Orientierungen für ein Lückenpaar [AouBa]. Diese Orientierungen 
eines Lückenpaares mit der assoziierten Aj auf dem +-G.P. {x, y, z} lauten: 


Tabelle 1 


| Ort der assoziierten | Ort der assoziierten 


Orientierung des Lückenpaares im Raum 4 
“ 
Orientierung 1 bzw. (1,0,0) ....... (z, y, 2) (x +1, y, 2) 
Orientierung 2 bzw. (0,1,0) ....... | (x, y, 2) (x, y + 1, 2) 
Orientierung 3 bzw. (0,0,1).. 2.2...) (x, %, 2) (x, y, 2 + 1) 
Orientierung 4 bzw. (—1,0,0) ...... | (x, Y, 2) (x —1,y,2) 
Orientierung 5 bzw. (0, —1,0) ...... (x, Y, 2) (x, y — 1, 2) 
Orientierung 6 bzw. (0,0,—1) ..... . | (x, 2) (x, y,2—1) 


Für die Ableitung der kinetischen Gleichungen sind maßgebend die bereits 
eingeführten Wahrscheinlichkeiten f ,(t, t), fz (t’, t), 9, (t, t). Weiter benötigt 
man noch die Wahrscheinlichkeit y (t, ¢) dafür, daß zur Zeit ¢ ein Lückenpaar 
7) Ein höheres Liickenaggregat kann bestehen aus drei oder vier (oder noch mehr) 
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iner 6 6 


Für die Wanderung der Lücken sowie für die Bildung der Lückenpaare sind 
nde | die entscheidenden Größen die Sprünge der benachbarten Gitterionen in 
hen | die Lücken hinein. Die verschiedenen Sprungarten sind in Tab. 2 zusammen- 
em gestellt: 


be- 
nd Tabelle 2 
en. 
4 Sorungzahl pro |__ Genaue Kennzeichnung des jeweiligen Sprunges 
~ | Zeiteinheit für | Abstand der | Abstand der 
ind | die jeweilige Zustand vor Zustandnach | “4 undB A.undB 
= y y o o o o 
lar. Lücke dem Sprung dem Sprung vor Sprung | nach Sprung 
t), a(V5) für A, A, ist isoliert A, ist assoziiert v5 1 ee 
be- }x(/3) für A, | A, ist isoliert A, ist assoziiert y3 1 Y ; 
i.e für A, A, ist assoziiert A, ist assoziiert 1 1 i a 
a(V5) fiir A, A, ist assoziiert | A, ist isoliert 1 V5 H 7 
a(/3) für A, | A, ist assoziiert | A, ist isoliert 1 v3 } : 
4% fir A, | A, ist isoliert A, ist isoliert 
ar- 6(V5) fir B, B, ist isoliert B, ist assoziiert V5 1 WER. u 
hl- g(V3) fir B, B, ist isoliert B, ist assoziiert y3 1 a sr 
A fiir B, B, ist assoziiert B, ist assoziiert 1 NW Pe 7 
(V5) fir B, | B, ist assoziiert | B, ist isoliert 1 : i 
(V3) für B, | B, ist assoziiert | B, ist isoliert 1 V3 of 
By fir B, B, ist isoliert | B, ist isoliert >y5 >YV5 ee 
m- 


In Tab. 2 haben die Sprungzahlen folgende Bedeutung: 


«(Y5) - At = Wahrscheinlichkeit dafür, daß während der Zeit At ein der Ap 

benachbartes Ag-Ion in diese Ag springt. Vor dem Sprung 

ten muß die A, isoliert sein und den Abstand /5 von der näch- 

1... sten By haben; nach dem Sprung muß die Ap an die BJ asso- 
Nr ziiert sein; also muß dann ein Lückenpaar vorliegen. 


Genauso sind auch die Bedeutungen der anderen Sprungzahlen & (y 3), gr 
usw. aus der Tabelle 2 zu ersehen. 


Sämtliche Sprungzahlen haben die Gestalt: &; = »; - exp| -; | mit: 
eits 
tigt &, = Schwellenenergie, welche beim Sprung A überwunden werden muß; 
mar k = Boltzmannkonstante; 
T = absolute Temperatur des Kristalls; 
v,ist proportional der Frequenz der Gitterschwingungen. 3 : 


| 


ann | mit der assoziierten A, auf dem +-G.P. r vorhanden ist. Es gilt offenbar: 
S80- 
i 
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Aus der formalen GI. (10) erhalten wir nun unter Verwendung der Sprung- 


zahlen in Tab. 2 die explizite Form für I y, (r, t) 


= -B(V5)- (@ + By t+ Leal) t+ fp (x4 
+ 2,y,2+ 1,t) + f, (2 + 2, y—1,z,t)+ fg(2+ 2,y,2-1,)) + 
+ fy (2, ys -B(V3)- (@ + + y+l,2- 
4+ + 
++ [ps (x + 1,y + 1,2, t) + po (x + 1, y—1, z, t) + yg +. 

+ 1,y,2+ 1,t) + yg (x + 1, yz—1,)) 
+f, (x +1, y,z,2)- 
+1,z—1)+ f,(z,y—1,2—1,8)] + B*- [ye (2, y, 2, t) + 

+ ps (x, y, 2, + ps (2, y, + (x, y, 2, t)] — y, (8, y, 2,8) - [4 4 
+4p*+4-a(/5) 


Die Ausdriicke fiir u Wo (t, t), a (r, ¢), ous (t, ¢) usw. erhält man genauso 


wie (11) aus der formalen Gl. (10). Übrigens kann man auch a Yo: 3 Ws USW. 
aus (11) durch Symmetrieperationen®) gewinnen. 

Die Gl. (11) sowie die 5 analogen Gln. für z Yo usw. sind nur gültig für das 
Innere (Volumen) eines Kristalls. Dagegen an Versetzungen, Mosaikgrenzen, 


inneren Oberflächen usw. müssen andere Gleichungen benutzt werden. Dieser 
„Oberflächenfall‘‘ wird hier nicht behandelt. 


é 
8) Differentialgleichung für nm”? Diffusionskonstante der Liickenpaare 


In den Ausdrücken (11) und ähnlichen für a Yo, a YW, usw. werden 

fy (to, t), fz (t + i, — ¢) in Taylorreihen®) entwickelt um den 
2 

Punkt r = {2, y,z} herum. Nach den quadrat. Gliedern ji werden die 

Reihen abgebrochen. Hierbei wird per definitionem gesetzt: f, (t,t) = 


8) So erhält man z. B. aus (11) den Ausdruck für 2 Wo (t, t), indem 1. x wird mit y ver- 


tauscht; 2. die Indizes von y, werden ersetzt durch (1, 0, 0) — 1,(0,1,0) -- 2, (0, 0,1) +3 
usw. Dann wird in (11) das Indextripel (1, 0, 0) vertauscht mit (0, 1, 0) usw. 

9) Die Funktionen f, (t, t), . . . sind zwar nur für isolierte Punkte, nämlich die Punkte 
tg.p. der Gitterplätze definiert. Jedoch kann man die Werte von f, (t,t),... für da- 
zwischenliegende Punkte r durch Interpolation definieren. Im Text wurde bei tq.p. der 
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6 6 
fs (t + Danach werden die y, durch f,, y = sowie deren 


Ableitungen dargestellt. SchlieBlich wird noch = y gebildet. Die Rechnungen 


sind recht mühsam, jedoch elementar und liefern als Ergebnis 


— 


(t,t) = 24 [0 (V5) + « (V3) + B(V5) + B(V3)) — 


(4, — = 8 [a (V8) + (V5) + +26 


In der Diff.-Gl. (12) ist = (V5) + (V3) + 
4 B(y5)) fı Isa [a(V 3) + a(y5 ) +6 +6 (V5)ly kennzeichnend 


für die chemische Reaktion: Ag + Ba @ [Ag Bo], wodurch Lückenpaare 
entstehen. 
ie Diffusi 
Der Term E Ch EER LS Ay (t, t) bedeutet die Diffusion 


der Lückenpaare [A-ı Bo]. Die Diffusionskonstante der Lückenpaare!®) 
muß offenbar sein: 


(13) 

u 2/ 3 + 


In kommt der Faktor ( 5) herein,weil im benutzten Koord.-System is ) 


Längeneinheit entsprach. In der Formel für D,.»., kann y im allgemeinen gegen a* 


werden den wegen (3; 5) kT 0,05 eVolt, Ae = 0,5 eVolt > 
a 


* 


= 4 (V3) + a(Y5) + 38(V3) + 5B 


grad f,. grad fy 


Der Term E vi stellt dar zusätzliche Glieder infolge räumlicher Inhomogeni- 
usatz 
täten von f,, fa, y;- Bei geringen Gitterkonzentrationen der Fehlstellen ist das Glied 
E ER klein. Es hat die Gestalt 
6 a* 
mas je a* + BY + y Ay) 4) Ay] - 


‚uso 
Sw. 
en, 
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len 
die 
te | 
la- 
ler 


= 
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Die Diffusionskonstanten D,, bzw. D,, der einfachen Lücken A, bzw. 

Bg betragen (s. Gl. (15)): Da, = 4%: (5) bzw. Dy, = 4 fo: Es gilt 
nun bestimmt wegen der geringeren Schwellenenergie x») < «*, By < ft. 
Daraus folgt aber: Die Lückenpaare [A By] müssen bedeutend beweglicher 
sein als die isolierten Lücken! 
0) Die Ditferentialgleichung für 


Bei der Berechnung von oh, wird als chemische Reaktion nur die Bildung von 


Lückenpaaren berücksichtigt. Dagegen andere mögliche Prozesse, wie z. B. die Erzeugung 
und Vernichtung isolierter Lücken an Versetzungen bzw. innere Oberflächen werden nicht 
behandelt. 
Unter dieser Voraussetzung erhält man mit f, (t,t) = NV P, (ty, t, t) für — f,(t, t) 

at 
die Beziehung"): 

6 
fy (t, 8) = (te, t) — + (to, t) a (te + 


at To j=1 t — 25) 


6 
— fa (to + +,—r)- O(r— 1,)%,2 _ 
G=1 (t— + iy) 
6 
7: fa (t,t) a he (tS, t) B (Ir, t|) : Or 19,2 (14) 
j=1 + iz) 


Aus der Gl. (14) wird nun der explizite Ausdruck für ote (t, t) abgeleitet. Ta ylor- 


reihenentwicklungen und weitere Umformungen liefern nun die Endformel: 


= 24 (V8) + (V5) + BYE) + 8) 4 
+4[a(V3) + a(V5) + 0(V3) + + + 


é 
15 
+ r 


Der Term [= RER in (15) hat eine ähnliche Gestalt wie der Term [& vl... in (12) 


Der Term E A = — 24 [a(Y3) + « (Y5) +B (V3) +B (VB) 


Reaktion 
fy (t,t) + 4[a (V3) + a(V5) +5(Y3) + 5(Y5)] y entspricht wieder dem Prozeß: 
A, + o Bo 
Das Glied = ha] = 409° Af, (t,t) stellt dar die Diffusion der A,-Lücken. 
Lat riffusion 
Die Diffusionskonstante der A,-Liicken ist hierbei: 


Eine zu (15) analoge Gleichung erhält man auch für Fr] Is (t, 2). 


1) Für f, (t,t + At) muß eine gegenüber (5) geänderte Beziehung benutzt werden. 
Die hinzukommenden Glieder tragen der chemischen Reaktion A, + B, = [As Bo] 
Rechnung. Man erhält diese Beziehung durch Betrachtung einer statistischen Gesamtheit 
gemäß (1). 
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«) Allgemeines über die Differentialgleichung für}, y(t,t), (t,t), Is) 


Aus den Diff.-Gln. (12), (15) sowie einer Diff.-Gl. für = fa (t, t) könnte 


man die Wahrscheinlichkeiten f4 (t, £), fz (t, £), y (t,t) berechnen, sofern die 
Anfangsbedingungen bekannt sind. Allerdings die praktische Lösung dieser 
Diff.-Gln. dürfte mathematisch recht schwierig sein, Außerdem ist das Ver- 


fahren nicht realistisch. Denn in den Diff.-Gln. für fa (rt, t) und 


“ fz (t, t) fehlen die anderen ,,chem. Vorgänge‘‘; insbesondere die Entstehung 


und Vernichtung der isolierten (einfachen) Lücken A, und BJ an Versetzungen, 
inneren Oberflächen usw. 


Allgemein kann man aber aus den Diff.-Gln. folgern: Die Wahrschein- 
lichkeiten f4 (t,t), fp, (t,t), p(t, ¢) haben das Bestreben im Verlaufe der 
Zeit t räumlich konstant zu werden, also unabhängig vom Ort r. Dies wird 
bewirkt durch die Diffusionsterme in den Diff.-Gln. (12), (15). 


Ferner kann man aussagen: Wenn zu irgendeinem Zeitpunkt t die Wahr- 
scheinlichkeiten f4 (t, r) = = y (t,t) = G 
räumlich konstant sind, dann sind die Wahrscheinlichkeiten f4, fg, y auch zu 
allen späteren Zeitpunkten t > t unabhängig vom Ort r, obgleich die chemische 
Reaktion Ag + Ba = [Ag Ba] noch ablaufen kann. Die letzte Eigenschaft 
kann bei anderen chemischen Prozessen fürdie einfachen Lücken ungültig werden. 
Dies ist der Fall, wenn nicht für alle Gitterplätze die Diff.-Gln. dieselbe Gestalt 
haben. 
d) Die kinetischen Gleichungen für die Bildung der Lückenpaare . | 
In dem Kristall des Volumens V seien an Fehlstellen enthalten: 


, also 


Nı4(t) = Anzahl der isolierten A, zur Zeit ¢ im Kristallgitter des Volu- 
mens V; 


Nz (t) = Anzahl der isolierten B,-Liicken zur Zeit t im Kristallgitter des Volu- 
mens V 


Nap (t)= Anzahl der Lückenpaare [Aj Bp] zur Zeit ¢ im Kristallvolumen V 


G = Anzahl der positiven Gitterplätze im Kristallvolumen V = Zahl 
der neg. G.P. 


Es gilt nun): 


(16) 
6 
2) In |» erfolgt die Summation über alle pos. G.P. des Kristallvolumens V. Die 
+ G.P. 
Relationen (16) können durch Betrachtung der statist. Ges. gemäß (1) abgeleitet werden. 
Bei hinreichend großem V ist Streuung vernachlässigbar. 


| 
T- 
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Durch Benutzung dieser Beziehung erhält man: a 
mogenit 
ai Nan = 24 [x(V5) + + + 24[0(V3) +] when 
+BV3))- ine +3 (Y3)+a(Y5)+ | avhingi 
+b (V3) +b (y5)] wirkung 
+ V5ex = | Isle + n| mit Aus 
1 3 Ionen « 
¢, = Vi (2%, + 4,),... usw. (17) weniger 
8 die Ag- 
= : +t, + 7,),... usw. einen K 
verschie 
Die </p> sind also mit Gewichtsfaktoren / , (r, t) behaftete Mittelwerte; deshalb 
müssen die (fs) irgendwie mit a zusammenhängen. Dies kann erreicht 
werden durch den Ansatz: Orient 
(t + | de, t)> = 5 (t) - (t) Is(t + t) = Net): 8 (1). 
(18) 
Somit erhält man als kinetische Reaktionsgleichung: Orienti¢ 
= = = 24 [(a(Y5) + B(Y5)) + («(V3) + B(V3)) ()- Orientii 
(19) 
Na Na _ 4[a(y3) + a(V5) + 6(V5)] - (0. Orienti 
Orienti 
Die zeitabhängigen Funktionen s, (¢), 8, (t) könnten aus den Gln. (12), (15) : 
berechnet werden. Daraus folgt: 0 < s, (t) < 1. Der Verlauf der Funktionen Orienti 
8, (t) hängt von den Anfangsbedingungen ab. Zu Beginn des Vorganges ist 
8, (0) = 0; im Verlauf der Zeit ¢ streben s,(t) > 1. Wenn zu einer Zeit r En 
rienti 
gilt f, (r.r) = a ‚so haben wir gemäß (17) s, (r) = 1. Aus den Gin. (12), 
(15, (17) folgt nun: Ist für einen Zeitpunkt t die Wahrscheinlichkeit f, (t,t) | Orienti 
unabhängig vom Ort, also f, (t,r) = u, so gilt für r und alle späteren 
Orienti 
Zeiten t stets: s,(¢) = 1 mitt >t. Im kinetischen Gleichgewicht z N ap =8 
ist sicher s, (¢) = 1. Deshalb erhält man als Massenwirkungsgleichung: 
D: 
a (V3) + (V3) + a (V5) + 5(Y5) m 
= 
6 [x (V3) + Ein di 
FEIERN 
Als Resumé kann man feststellen: Die atomistische Behandlung der Be 
Bildung von Lückenpaaren liefert für die kinetische Reaktionsgleichung 13) 
eine gewisse Modifikation (19). Diese ist bedingt durch die räumliche Inho- = 
>) 


ite 


er 
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mogenität der Wahrscheinlichkeit f, (t,t). Dadurch kommen in die kineti- 
schen Konstanten (19) noch die zeitabhängigen Funktionen s, (¢), s, (¢) herein: 
0 < s, (t) < 1. Die kinetischen Konstanten in (19) zeigen also noch eine Zeit- 
abhängigkeit. Gewonnen wurde weiterhin der Ausdruck für die Massen- 


wirkungskonstante (20). 


IV. Bildung von Assoziaten [AgoSg] in Silberhalogeniden 


Aus isolierten Silberionen auf Zwischengitterplatz Ag, und isolierten Sg - 
Ionen entstehen die Komplexe [Ag.,Sg]. Die Sg -Ionen sind auf einigen 
wenigen —-G.P. mit Ortsvektor r, = {x,, y,, 2,} fest eingebaut. Dagegen 
die Ags -Zwischengitterionen sind frei beweglich. 

Bei einem Komplex [Ag¢6S¢] sitzt ein AgS-Ion in unmittelbarer Nachbar- 
schaft vom Sg -Ion auf einem Zwischengitterplatz (Z.G.P.). Es gibt nun für 
einen Komplex [Ag>S¢] mit dem Sg -Ion auf —-G.P. 1, = {2,, y,, z,} acht 
verschiedene Orientierungen im Raume: 


Tabelle 3 


( SSO- *S ASSO- 
Orientierung des Assoziates [Ag;S/,] im des Ort des 
gis ziierten zllerten 
Raume 
Ag;-Ions Sj,-Ions 

ae 1 1 1 
RER 1 1 1 
1 1 1 
Orientierung (+,—,+)oder 3..... Wis y 
DZ 1 1 1 
Orientierung (+, —,—) oder . . Len 2, 
1 1 1 
Orientierung (—, +, +) oder 5..... 2. 
Es 1 1 1 
Orientierung (—,+,—) oder 6... . . | u. - 
2 2 2 

1 1 1 
Orientierung (—, —, +) oder 7. .... | 
1 1 1 
Orientierung (—, —,—) oder 8. . . .. 2:08, %— 


a) Die Bewegung der Ag;-Ionen auf Zwischengitterplatz (Z.G.P.) 


Die Fortbewegung der Ag.;-Ionen beruht wieder auf Sprüngen. Nach 
J. Teltow!?-15) ist der Sprungmechanismus der Ago-Ionen folgender: 
Ein dem Ag; -Ion benachbartes Ag@-Ion auf G.P. springt auf einen umgebenden 
leeren Z.G.P., und das ursprüngliche Ag;-Ion springt auf den leer gewordenen 


) J. Teltow, Ann. Physik 5, 63 baw.71(1949). 6 
“) Seitz, Acta Cryst. 8, 355 (1950). mit: 
») J.E. Hove, Physic. Rev. 102,915 (1956). 
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zliert sein, mit Abstand + 3 vom Sg -Ion. 
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+-G.P. Im Endzustand solch eines gekoppelten Sprunges sitzt das ursprüng- 
liche Ag,-Ionaufeinem +-G.P. und ein anderes Ag- Ion auf dem neuen Z.G.P. 
Außerdem vermag das Ag¢-Ion auch direkt auf einen benachbarten Z.G.P. 
zu springen. 

Gemäß diesem Mechanismus gibt es drei verschiedene Sprungarten A, B,C 
mit verschiedenen Sprunglängen V1, V2, V3. Als charakteristische Größe 
wird eingeführt der Sprungvektor 3 gemäß: Sprungvektor $ = Ortsvektor vom 
Ago nach Sprung — Ortsvektor vom Ags vor Sprung. Die A- Sprünge mit 
Sprunglängen [84] = -y1 haben die Sprungvektoren: $4 =i, 
B43 = t,, 844 = — ty» 845 = — ip 84g = —?,. Für die B-Sprünge mit Sprung- 
lange = Y2 lauten die Sp, = t, + ty Spo = i, +i, 
Sp, = ti, + 7, usw. Es gibt insgesamt 12 Sprungvektoren # 8p. Die anderen 
entstehen durch Vorzeichenvertauschung aus $,,, Sp, 83; Für die C- 
Sprünge mit Sprunglinge |3c| = //3 haben wir als Sprungvektoren $e: 
entweder im Falle 1 


oder im Falle 2 
> . . . . . . . . . 
+i,—i. 
u z 


Diese Sprungarten sind also maßgebend für die Bewegung der Ag;-Ionen. 
Jedoch gilt dies nicht!%)1?) für radioaktive Isotope als Ag&-Ionen. Denn 
diese sind markiert, und man muß stets ein und dasselbe Teilchen betrachten. 
Während hier verschiedene Ag*-Ionen an einem Sprung beteiligt sind und es 
nicht auf ein bestimmtes Agt-Ion ankommt. Be 


b) Die Differentialgleichung für die Assoziate [Ag;S/,] ‘ 
Für die nun folgenden Rechnungen werden als maßgebende Größen be- 
nutzt: 
fo (t, t) = Wahrscheinlichkeit dafür, daß zur Zeit t auf dem Zwischengitter- 
platz (Z.G.P.) mit Ortsvektor r ein Silberion Ag, sich befindet; 
f,(t,.¢) = Wahrscheinlichkeit dafür, daß das auf dem —-G.P. r, eingebaute 
Sg -Ion zur Zeit ¢ isoliert ist; 
y,(t,. t) = Wahrscheinlichkeit dafür, daß das auf dem —-G.P. r, eingebaute 
Sg -Ion zur Zeit ¢ mit einem Ag,,-Ion gemeinsam ein Assoziat 
[Ag5S¢] der Orientierung j bildet. 
Die bei der chemischen Reaktion Ag 5 + SZ [Agö5S&] in Frage kommenden 
Sprungmöglichkeiten sind in der nachstehenden Tab. 4 enthalten. 
In dieser Tab. 4 bedeuten wiederum die Sprungfrequenzen z. B.: 
ßz: At = Wahrscheinlichkeit dafür, daß im Zeitelement At ein Ag;-Ion einen 
gekoppelten Sprung der Sorte B ausführt; hierbei mußte vor dem 


Sprung das Ag,-Ion isoliert sein und den Abstand 2 yıı vom 


-Ion haben; nach Sprung muß das neue Ag,-Ion asso- 


16) J. Bardeen u. C. Herring, in nearly Cope, S. 281, 
Herausgeber: W. Shockley u. a. 
17) Robert J. Friauf, Physic. Rev. 105, 843 (1957). en we oT, ae 
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Tabelle 4 
Genaue Kennzeichnung des jeweiligen Sprunges 
Sprungfrequenzen für Zustand des | Zustand des aa _ _ 
Ag;-Ionen auf Z.G.P. Ag:-Ionsvor Ag:-Ionsnach | ‚?w’schen zwischen 
Ag; und Si, | Ag; und Sj 
prung nn. vor Sprung nach Sprung 
| | 
27 
«, für Sprunglänge yı isoliert | isoliert >F -y27 
a, für Sprunglänge 2 isoliert isoliert S y2? | > V2 
% für Sprunglänge ) 3 isoliert | isoliert > + y27 > y27 
| Am 
B, für A-Sprung isoliert _ assoziiert | > yi1 3 “V3 
Ba für A-Sprung assoziiert | assoziiert | —'V3 V3 
1 / 1 
B, für B-Sprung isoliert assoziiert > yı 5 -V3 
| 1 
für B-Sprung isoliert assoziiert -y1ı9 + -V3 
Bs für B-Sprung assoziiert assoziiert v3 | 
8, für C-Sprung isoliert assoziiert “Vil + V3 
en 1 
ße für C-Sprung isoliert assoziiert | °V3 
b, fiir A-Sprung assoziiert isoliert yıı 
A pee 1 
b, für B-Sprung assoziiert | isoliert 5 -/3 + “Vil 
= 
b, für B-Sprung assoziiert isoliert 5 -Y3 5 -y19 
b, für C-Sprung assoziiert isoliert : v3 
b, für C-Sprung assoziiert isoliert .y27 
_ fA 
Die Sprungzahlen pro Z.E. haben die Gestalt: a, = 2-¥y-e #+ 


e4 
-ım 


= va: exp(- x = exp (- usw. 


Mit den eingeführten Wahrscheinlichkeiten fy, f, y; und Sprungfrequenzen 


erhält man die Beziehung?) 


18) Gl. (21) wird wieder abgeleitet mittels der statist. Gesamtheit gemäß (1). Die 


Einheitsvektoren ¢, haben die Gestalt: 
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| herum 


Aus dies 
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Aus dieser formalen GI. (21) erhält man nun für £ Yır.ı 
pliziten Ausdruck: 
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Die Ausdriicke fiir a Pett) (tt) usw. können aus Gl. (22) erhalten 


werden durch Symmetrieoperationen oder direkt aus der formalen Gl. (21). 
e) Die Differentialgleichungen für die Wahrscheinlichkeiten fy (t. £) der Ag;-Ionen 


Für einen Zwischengitterplatz (Z.G.P.), in dessen Umgebung kein Sg - 
Ion auf einen en P. eingebaut ist, gilt die Diff.-Gl. 


ot ( )= fo ( e -TB,C u (rt to) Su 
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Aus dieser formalen Beziehung wird zuerst die explizite Darstellung abgeleitet. 
Anschließend erfolgt eine Taylorreihenentwicklung von f, (t,, £) um den Ort 
therum. So bekommt man schließlich die Formel: 


(t, t) = [l-a, + 4a, + 20.]: Ah (t, 2). (23) 


Gl. (23) gilt nicht für radioaktive Isotope; denn dann ist das einzelne Ag’-Ion 
markiert”). 
Gl. (23) gilt für Z.G.P. r, in deren Umgebung kein Störion Sg eingebaut 


ist. Die Diff.-Gl. für ale wird dagegen modifiziert, wenn in der Nachbar- 


schaft des betreffenden Z.G.P. ein Sg -Ion a ist auf —-G.P. Beispiels- 


weise erhält man so für einen Z.G.P. + + + +} die Diff.-Gl.: 


+ be Yo lty + (tet) + 

+ (tes + 94 4,4 (tert) — [Bg + 2 + Bel: fs (ty, - 


(24) 


+ 


to] we 


1 
Die Werte fy, (x He t) usw. wurden per definitionem ein- 


geführt und zwar durch Interpolation der umgebenden Werte. Auf solch einen 
1.6.P. + ty 3 es) ist ja das Ag~-Ion assoziiert, deshalb hat dort f, + 


1 P 
3/3 en) keine physikalische Bedeutung. - 
Analoge Diff.-Gln. wie (24) erhält man auch für 


fa (2, + Ys t) ’ Ys +% +2 t) usw. 
Die Gln. (23), (24) zeigen also: Nicht auf sämtlichen Z.G.P. gelten die- 

slben Diff.-Gln. für die Ag,-Ionen. Dies steht im Gegensatz zu den Gln. für 

die A4- und By-Liicken (s. Gl. (15)). 
Aus den Diff.-Gln. (23) und (24) folgt nun: Die Wahrscheinlichkeiten 

h(t, t) werden erst dann unabhängig vom Ort rt, wenn die chemische Reaktion 

Ags + Sq = [Ag6S¢] ihr kinetisches Gleichgewicht eingestellt hat. Auch dies 

steht im Gegensatz zur Lückenpaarbildung. 
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I 
ta d) Kinetische Gleichung für die Bildung der Assoziate [Ag; S3]. h: der 
Kinetische Konstanten. Massenwirkungskonstante 
ie Benötigt wird nun die Anzahl der jeweiligen Störstellen: " iz 
N,(t) = Anzahl der isolierten Ag¢-Ionen auf Z.G.P. im Kristallgitter def. 
Volumens V zur Zeit t; Die kin 
N,(t) = Anzahl der isolierten Sg -Ionen im Kristallvolumen V zur Zeit t, 
Nass (t) = Anzahl der assoziierten Sg -Ionen im Kristallvolumen V zur f 
Zeit t; also die Zahl der Komplexe [Ag¢S@]. Nass (1 
Es gelten nun die Relationen: 
t) = s (t,, t). he dui 
eiten 
Die (22) nun in (29) eingesetzt. Dadurch erhält man it durc 
längerer aber elementarer Rechnung die Formel: Im kin 
d 
qi Nas (t) = 24(B, + + Be)‘ (fo zyu t)). N, (t) + psetz. 1 


+ (ho (x. +5V19- ex, t) -N,(t) + 8B, + 


27 t)) N, (t)— (3b, + 6b, + 36, + 365+ do]: 


Hierbei werden die Mittelwerte (f,> gebildet mit den f, (r,, ¢) als Gewichts m ist 
en 


faktoren. Daher gilt: zen € 
1 1 24 1 
fo (ts + zyu "Co, t)) = = fs (tg, 8) - fo (x, +2 Yıil-e, t) usw.| ¢) Fo 
y = In Ag! 
Die Einheitsvektoren e,, e* und e, haben die Gestalt: + Zusat 
= (34, +4, +4,), = = („+ 32, + 7,),...es gibt e, bis e, us ex 
ub. Die 


1 


-(3i, +4, + 3%,),...es gibt e* bis es itterkon 


1 
yı9 


ya &= („+i,—i,), ...es gibt e, bis 
eführt 
Die Mittelwerte </,> müssen irgendwie mit =. zusammenhängen. 


_ Dem wird RUE getragen durch den Ansatz: 


/ 0 h wi 

+ 2 zVi9- en ))= (fo (x + “Nol t auch 
Hierbei: Hert da: 


Z = Anzahl der Zwischengitterplätze im Kristallvolumen V; Z = 2G. 
Die zeitabhängigen Funktionen s, (¢) erfüllen die Bedingung: 0 < s, (t) <1 
= zeitlicher Verlauf hängt sehr von den Anfangsbedingungen ab; d. h. vom 


®) Chi 


i 
| 
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itt, der Sg -Ionen sowie der Anfangsv erteilung der Ago-Ionen. Zu Beginn SHE 
is Vorganges ist s, (0) = 0 und im Laufe der Zeit geht s, (t) > 1. Der Wert _ aah 
= 1 ist exakt erreicht, wenn das kinetische Gleichgewicht der Reaktion Bs! bs 

+ Sg = [AgoS8¢@] eingestellt ist. Dies steht im Gegensatz zum Ver- 
tabi der s, (£) bei der Lückenpaarbildung. 


Die kinetische Reaktionsgleichung lautet somit: 


eit t; 
zur 


EN (t) = [24- 8, (t) (By + + Bo) + 24> 8 (t) - Bg 83 (8) - Be] 


BL [3b, +65, + 3b, + Nass (). 


(25)] Die kinetischen Konstanten zeigen also eine zusätzliche Zeitabhängigkeit, a 
Iche durch die räumlich inhomogene Verteilung der Aufenthaltswahrschein- | 
hkeiten bedingt wird. Weiterhin wurden die kinetischen Konstanten ex-  _ 
nach fizit durch die Sprungfrequenzen ausgedrückt. 


Im kinetischen Gleichgewicht: = N ass (t) = 0 gilt das Massen wirkungs- 


rsetz. Dafür ist bestimmt s, — 1 erreicht. 


a 
N, N, _ + 6 by + + 3 be + be a7) 


Z 8[3 6, + 6B, + 3B, + 3 Be + Be) 


38 (t). 
=” ist auch die Massenwirkungskonstante in (27) durch die Sprungfre- 


nzen explizit dargestellt worden. 


usw.] e) Folgerungen aus dem Massenwirkungsgesetz fiir die Assoziatbildung 


In AgBr mit etwa 0,01 Mol-% Ag,S-Zusatz zeigt die Lichtabsorption !?) = = 
t Zusatzkristalle bei —190 °C eine Bande bei 493 mu. Diese Bande wird | 
f isolierte -Sg-Ionen zurückgeführt. Nach der Dispersionsformel tet 
2 us experimentell, daß mehr als 1%, der Sg-Ionen bei —190 °C isoliert sein <= 
ub. Dies soll nun mit dem M.W.G. (27) konfrontiert werden: Werden die 


e% piterkonzentrationen 
N, N, N 
Pg: @’ » CAss = mit Z=2G 
peeführt, so wird aus (27) die Beziehung: wl 
yen. 
"Cy 2 136,+66,4+36,+36,+56 
CAss 8 [364 +6 Be +3 + 3 Bo + Be 
Die Schwellenenergien in den Sprungfrequenzen sind zwar unbekannt, ae 
bch wird bestimmt gelten: -“ we FP mit Ae > 0,30 eVolt. Analoges 
: hi auch für die anderen Sprungzahlen. Die Kristalltemperatur 7’ = 80 °K ae: =: 
Hert dann 44, also e# 7,8- 10%, Somit folgt aus dem M.W.G 
von ) Chr. Volk, Z. Physik 138, 623 (1954). aeg EN 


= 
(26) 
| | 


NER, Annalen der Physik. 7. Folge. Band 5. 1960 


wenn die anderen Sprungzahlen als klein vernachlässigt werden: 


— 18 19-20 
< 48. 10-, 


CAss 


krepanz beseitigt!! 


Eichwalde bei Berlin, Stubenrauchstr. 90. 


y 
Bei der Redaktion eingegangen am 14. Mai 1959. 


: * 


Nun ist cass © Czusatz = 10%. Wird noch c, = c, angenommen, dann ergi 
sich schließlich c, < 10-*. Dies steht im krassen Gegensatz zum experim 
tellen Befund von c, > 10%. Das Massenwirkungsgesetz (27) kann 
schwerlich modifizieren. Jedoch kénnen noch andere ,,chem. Reaktione 
a der Agü-Ionen stattfinden. Beispielsweise könnten die Ags-Ionen an be 
A Er deren Stellen wie Versetzungen, Mosaikgrenzen usw. im Kristall bevorz 

angelagert sein. Falls solche Stellen als Haftstellen wirken, wäre die Di 
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Die Richtungsabhängigkeit der Magnetostriktion 
tio’ 


Von W. Döring und G. Simon 
\ 


N 
In Einkristallen läßt sich die Magnetostriktion durch einen sy mont by 
weistufigen Tensor beschreiben, dessen Komponenten von ‘der Magneti- 
erungsrichtung abhängen. Diese Abhängigkeit wird mit Hilfe eines gruppen- 
eoretischen Verfahrens für alle Kristallklassen in eine solche Form gebracht, 
ß sie der Kristallsymmetrie genügt. Die Resultate sind nicht nur auf die en 
agnetostriktion, sondern auch z. B. auf den Widerstand in magnetischen 
inkristallen und in der Theorie der Gitterschwingungen anwendbar. 


I. Problemstellung 


In einer früheren Arbeit!) wurde die allgemeinste Gestalt der Kristall- er be 

ergie eines ferromagnetischen Materials fiir einen Kristall beliebiger Symme- 
berechnet. Hier soll dieselbe Aufgabe für die Magnetostriktion gelöst 

rden. Die Verzerrung des Gitters durch die spontane Magnetostriktion 

ird durch einen symmetrischen Tensor beschrieben, dessen Komponenten 

HW. (i, k = 1, 2, 3) Funktionen der Richtungskosinusse x, der Magnetisierungs- 

iehtung sind. Die relative Dehnung A = (Al/l),, einer Strecke der Länge / 

Richtung mit den Richtungskosinussen ß, ist dann durch den Ausdruck 


Jede Symmetrieoperation des Kristalles führt eine beliebige Richtung x; 
neine andere x; über, was sich durch eine Beziehung der Gestalt 


3 
sdrücken läßt. Darin ist R,, eine reell orthogonale Matrix. Die Symmetrie 
Kristalles besagt für die Magnetostriktion, daß die in Richtung Bi gemessene 
gnetostriktive Dehnung fiir eine in Richtung x; dieselbe ist 
vie für das Richtungspaar «,, ß, 
A (oxi, Bi) = By), 
obei x; und ß; durch dieselbe Operation (2) aus x, und ß, hervorgegangen _ 
ind. 
Es soll nun die allgemeinste Gestalt der Matrixkomponenten 4,, (x,) 
kstgestellt werden, welche mit dieser Symmetrieforderung verträglich ist. 


I) W. Döring, Ann. Physik (7) 1, 102 (1958) (im folgenden als I zitiert). 


1 


= 


DESO q 
vorzu 
ie Di 
| 
| 
| 
| 
| 
| 
| 
2 
> 
| 
| u 
| 
| 
| 
| ‘ 
| 
1 


alle Kristallklassen gegeben werden. 
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Fir einen kubisch-holoedrischen Kristall wurde die Lésung bereite von Becks 
und Döring?) angegeben. Für hexagonale Kristalle hat Mason?) diesel 
Aufgabe gelöst, jedoch unter der einschränkenden Annahme, daß die 4 
nur Polynome der «; von höchstens 6. Grade seien. Das Verfahren y 


Mason läßt sich nur schwer auf kompliziertere Funktionen übertragen. | 


folgenden wird für die Magnetostriktion durch eine ähnliche gruppentheor 
tische Methode wie früher für die Kristallenergie die allgemeine Lösung fi 


Il. Die Symmetrieforderungen für die Komponenten des Verzerrungstenson 
Um die Symmetrieeigenschaften der A,, zu ermitteln, formen wir d 
Ausdruck (1) etwas um. Aus den 6 linear unabhängigen Richtungsfunktio 
P;ß, bilden wir 6 orthogonale (und normierte) Linearkombinationen 4, (f 
v—1,2,..,6. Diese neuen Größen sind Funktionen der Richtung ;. Unte 
dem skalaren Produkt (f, g) irgend zweier Funktionen f und g der u; 
P; soll das Integral über die Einheitskugel verstanden werden: 
= ff f* gdQ (dQ2 = Flächenelement). (4 
Einheitskugel 
Jede Symmetrieoperation R definiert nun einen auf jede Funktion f der Rich 
tung anwendbaren Operator P, durch die Definitionsgleichung®): 
Pr f (Bi) = f (Bi)- (i 
Bei der Anwendung von P, auf eine der Funktionen g, muß eine Linear 
kombination der g, entstehen, denn jede Drehung, Spiegelung oder Dreh 
spiegelung verwandelt eine homogene quadratische Funktion der ß, in ein 
gleichartige Funktion. Also gilt: 


= & Dur{R) | 


re 

Die Matrizen D,, (2) bilden eine Darstellung der Gruppe der Symmetrie 
operationen. Es soll hier angenommen werden, daß durch passende Wahl de 
gp, diese Darstellung bereits in ausreduzierter Form vorliegt, doch ist das ; 
das Folgende zunächst nicht erheblich. Bei den späteren Rechnungen ist e 
deshalb von Vorteil, weil man die irreduziblen Darstellungen aller kristalle 
graphischen Gruppen kennt. 

Die Matrizen D,, sind unitär. Jede Symmetrieoperation eines Kristall 
führt ein Flächenelement der Einheitskugel in ein gleich großes anders 
Flächenelement über; deshalb läßt der Operator P, das skalare Produkt 
(f, g) invariant: 


(Pref, (rg) (i 
Wendet man diese Beziehung auf 9, und y, an, setzt (6) ein und benutzt di 
> DR, (R) Dus (R) = (8) 
x=1 


kur ie 
*) R. Becker u. W. Döring, Ferromagnetismus 1939, S. 274. 
3) W.P. Mason, Phys. Rev. 96, 302 (1954). i 

4) Auch die Definition P, f(ß}) = f{ß,) wird manchmal verwendet. ‘if. 
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Durch Verwendung der Funktionen 9,(ß,) an Stelle der Produkte 8; 8, 
geht (1) über in 


A (a; Bi) = = 9% (B;). (1’) 


Die darin auftretenden Piinketenei yr (&,) sind gewisse Kombinationen 
der A;,, welche unter Umständen komplex sind, nämlich dann, wenn die 9, 
es sind. 

Für die zu den Koeffizienten y} konjugiert komplexen Koeffizienten y, 
gilt folgender Satz: 

Die y, spannen dieselbe Darstellung der Symmetriegruppe auf u die p,. 


Beweis: Man setze in (1’) die Definition (5) ein: i pa 
Bi) = P, vr Pp Pu (B,)- 
Benutzung von (6) ergibt: x 
6 wd 
40, Bj) = P, vr (,) (R) - Po (B,). 


Wenn man nun die Beziehung (3) ausnutzt und (1’) einsetzt, erhalt man: 
6 


6 
= 


u 


Diese Gleichung multipliziere man mit x D* p*, integriere ß; über die Einheitskugel 


und nutze die Beziehungen (8’) und (8) aus; das ergibt: 


6 
(&) = ~ Dr vr 


Da P, ein reeller Operator ist, gilt diese Beziehung auch für die zu y* konjugiert-komplexe 


Funktion y,: 
= = D,,(R) y,(&), 


womit die obige Behauptung bewiesen wäre. 


Die Aufgabe, die allgemeine Gestalt der Magnetostriktionsfunktion zu 
finden, ist damit auf folgende Weise zu lösen: Man bildet die Linearkom- 
binationen 9, der ß,ß, so, daß sich die sechsdimensionale Darstellung D,, 
in ausreduzierter Form aus irreduziblen Darstellungen zusammensetzt. Man 
sucht dann zu jedem darin enthaltenen irreduziblen Bestandteil die allgemeine 
Form des dazu gehörigen Funktionenkomplexes y, (&,). Setzt man diesen in 
(1’) ein, so hat man die allgemeine Gestalt von A («,, ß,) gefunden. Drückt 
man darin die gy, wieder durch die ß, 8, aus, so ergibt sich auch die allgemeine 
Form von jeder mit der Kristallsymmetrie verträglichen Tensorkomponente 
4 


ik 


III. Die allgemeine Form der Funktionen zu einer irreduziblen Darstellung 


Im folgenden seien D,,(R) (u,v = 1,2,..,p) die Komponenten der 
unitären Matrizen einer p-dimensionalen irreduziblen Darstellung einer der 
32 Kristallgruppen. Wir suchen die allgemeine Gestalt der Funktionen y,, (&,). 
welche (6) befriedigen. Dabei geniigt es, eine dieser Funktionen zu ermitteln, 
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etwa y,, denn bei bekanntem y, folgen die übrigen daraus eindeutig wegen?) 
p 

Yu Din (R)- Pry (9) 

(g = Ordnung der Gruppe). Die Bestimmung der allgemeinen Form von 

kann nach zwei verschiedenen Methoden ausgefiihrt werden. Die erste macht 

Gebrauch von den Projektionsoperatoren auf eine bestimmte Zeile einer Dar- 

stellung. Ist F (x,) eine beliebige Funktion der Richtung, so ist 


= & Di, (R) Pp F (&,) (10) 


für beliebiges » eine Funktion zur 1. Zeile der betrachteten Darstellung. Nach 
Melvin®) kann man diese Formel für spezielle Darstellungen erheblich ver- 
einfachen. Jede reguläre Funktion der Richtung kann man nach Kugel- 
funktionen 

Yim =( 42 |m))! d cos (11) 


entwickeln. Setzt man Y, „ an Stelle von F in (10) ein und läßt y, 1, m alle 
möglichen Werte durchlaufen, so erhält man im Prinzip unendlich viele Funk- 
tionen /, zur ersten Zeile der betrachteten Darstellung. In den meisten Fällen 
ergibt sich jedoch null. Leider sind die von null verschiedenen Funktionen 
nicht linear unabhängig. Wählt man eine Folge von linear unabhängigen 
Funktionen aus, so kann man die allgemeinste Funktion als unendliche Reihe 
nach diesen Funktionen mit willkürlichen Koeffizienten schreiben. 

Für die praktische Anwendung solcher Reihen ist es bequem, wenn die 
Entwicklungsfunktionen orthogonal und normiert sind. Im allgemeinen ist 
aber das Orthogonalisieren und Normieren der linear unabhängigen Funktionen 
ziemlich schwierig. Einfach ist es nur bei denjenigen Kristallgruppen, bei denen 
alle Symmetrieoperationen aus Drehungen um die z-Achse, der Spiegelung an 
der x—y-Ebene und y—z-Ebene und der Inversion am Koordinatenursprung 
zusammengesetzt werden können. Das ist bei allen Kristallgruppen mit Aus- 
nahme der 5 kubischen 7, T,,O, T,und O, der Fall. Diese genannten Symme- 
trieoperationen reproduzieren nämlich alle Y,„ bis auf einen konstanten 
Faktor oder verwandeln Y, „ in Y,_„- Je nach der Kristallgruppe und der 
betrachteten Darstellung sind dann die von null verschiedenen Funktionen, 
die man beim Einsetzen der Y, „ in (10) erhält, abgesehen vom Zahlenfaktor, 
wieder die Funktionen Y, „ selbst oder die Funktionen ones 
C.0= Vio 


3 
und = 4 


= — Y,,m) 


In jedem Fall läßt sich leicht berechnen, fiir welche / wie m ‚eich null oder 
diese Funktionen ergeben. Ein Orthogonalisieren und Normieren erübrigt 


5) E. Wigner, Gruppentheorie und ihre Anwendung auf die Quantenmechanik der 
Atomspektren. Braunschweig 1931, S. 120. 
& 8) M. A. Melvin, Rev. Mod. Physics 28, 18 (1956). I NR 
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sich dann, denn diese Funktionen sind normiert. Ferner sind sowohl die Y, „ 
untereinander als auch die Funktionen C, „ und 8, „ untereinander ortho- 
gonal. 


Für die kubischen Gruppen wurde ein anderes Verfahren zur Ermittlung 
der allgemeinen Form der Funktionen benutzt, welches in einer Verallge- 
meinerung des in (I) benutzten Verfahrens zur Bestimmung der allgemeinen 
Funktionen zur identischen Darstellung besteht. Dieses Verfahren verzichtet 
auf Normierung und Orthogonalisierung der Entwicklungsfunktionen und er- 
strebt statt dessen die Herstellung einer möglichst einfachen, geschlossenen | 
Form der gesuchten Funktionen. Hierbei ermittelt man zunächst, wieviel vr 
linear unabhängige Linearkombinationen von Kugelfunktionen /-ter Ordnung 
es zu einer Zeile der betrachteten Darstellung gibt. Diese Zahl ist gleich der 
Häufigkeit des Vorkommens dieser Darstellung in der (22 + 1)-dimensionalen 
Darstellung, welche diese Kugelfunktionen von der Kristallgruppe aufspannen. 
Diese Anzahl kann man nach der Methode von Bethe leicht aus den Charak- j 
teren der Darstellungen ermitteln. In Tab. 1 ist das Resultat für die zwei- oe 
dimensionale Darstellung E der Gruppe O angegeben. Die zugehörigen Kom- 
binationen von Kugelfunktionen /-ter Ordnung kann man als homogene 
Polynome /-ten Grades der a, schreiben. Verschiedene solche Polynome 
ten Grades zur gleichen Zeile der betrachteten Darstellung sind linear unab- 
hängig. Außerdem sind sie linear unabhängig von allen Polynomen gleicher 
Art mit anderem Grade l’ +1; zu diesen sind sie sogar orthogonal. Wenn 
man auf Orthogonalität verzichtet und Bildung von irgendwelchen linear unab- 
hängigen Linearkombinationen als unwesentliche Abänderung ansieht, kann 
man sie also auch charakterisieren als linear unabhängige Polynome /-ten 
Grades der x, zu einer Zeile der betrachteten Darstellung, welche nicht aus 
gleichartigen Polynomen (/ — 2r)-ten Grades durch Multiplikation mit 1 = 


(af + «3 + a3)" und Linearkombination gebildet werden können. Wir suchen er 
mn solche Polynome zu bilden. Hat man 
wviele gefunden, wie die obige gruppentheo- Tabelle 1 


retische Abzählung angibt, so ist klar, daß alle Die Zahl der linear unab- 


weiteren von den gefundenen linear abhängig hängigen, homogenen Poly 
sind nome vom Grad / zur zwei- 


dimensionalen Darstellung 
An dem genannten Beispiel der Darstellung E der Kristallgruppe O 


von O sei das Verfahren weiter erläutert. Dabei 


oll im folgenden immer gleich ein Paar von zu- I= 12k 2k 4 
für das Funktionenpaar zu / = 2 findet man 12k +3 ok 
durch Anwendung von (9) und (10) oder auch ein- 1=12k+4 2k 4 is 
facher durch Raten das Funktionenpaar I=12k+5 2k+ 1 DEN 
12k +6 2k +1 
2 2 l=12k +7 2k+1 
fey = — 45, = 7, (at + 08 — 208) 
1=12k+9 2k +1 
und ebenso zu / = 4 1 = 12k 4 10 Ok 
= of — of, haa = 7, (ott a$ — 2 af). 12k +11 2k+2 
€ 
Aus jedem dieser Funktionenpaare erhält man ein weiteres, wenn man = 
ide Funktionen mit derselben, gegenüber allen Symmetrieoperationen == 
invarianten Funktion multipliziert. Die einfachsten, invarianten Funktionen =| 
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und p = at ab 


= 01 — 03) (a3 — 03) (03 — of). 


4 Eine beliebige invariante Funktion läßt sich nach (I) mit 2 willkürlichen 

Funktionen von s und pin der Form F (s, p) + q@ (s, p) schreiben. Multipli- 

ziert man die obigen Funktionenpaare mit s’ pi und gs p (i, 7 = 1, 2,...), 80 

2 man unendlich viele homogene Polynome zur betrachteten Darstellung, 

_ aber nicht alle, nämlich z. B. keine zu 1 = 5 und / = 7. Für diese Potenzen 
Bia man leicht unmittelbar 


1 P ‘ 
= Vs + — 208), = 04 (03 — af) 


l ° 
hex = + — 2 06§), foo = 05 — arf). 


Erweitert man diese Funktionenpaare ebenfalls mit s‘ p und gs’ p’, so 
erhält man insgesamt mehr homogene Polynome zu jeder Zeile der betrachteten 
Darstellung, als die Zahl der linear unabhängigen beträgt. Zwischen ihnen 
müssen also lineare Abhängigkeiten bestehen. Mit einiger Mühe findet, man 
diese heraus. Sie lauten: 


P) 


9 


Diese Gleichungen besagen: Alle Funktionenpaare, die man durch Multi- 
plikation der 4 Funktionenpaare zu / — 2, 4, 5 und 7 mit q - sö - pP (i, 7 beliebig, 
ganz) erhält, sind linear abhängig von denen, die man bei Multiplikation der- 
selben Funktionenpaare mit allen s* p” (k, n beliebig) erhält. Die zuletzt ge- 

nannten Funktionen ohne den Faktor q sind alle linear unabhängig. Die Hälfte 

im nämlich in x, gerade, die andere Hälfte ungerade. Drückt man ferner in 
ihnen a, durch «, und a, aus, so lautet das Glied mit der höchsten Potenz 
von a, innerhalb einer dieser Hälften bei allen verschieden. Ordnet man sie 
nach der Summe ! der Potenzen von a, a, und «3, so findet man für jedes | 
Übereinstimmung ihrer Anzahl mit der oben gruppentheoretisch ermittelten 

Zahl. Also sind damit alle linear unabhängigen Polynome gefunden. Ein 
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durch Multiplikation jedes der gefundenen Paare mit einem Zahlenkoeffi- 
zienten und Summation. Dieses hat daher die rasan 
fh = (af — 08) F, (8, p) + (at — af) F, (8, p) + % Xy (AT + — 208) 
Fs (8, p) + 04 (xf + af — 209) F, (8, p) 
lichen} 7s 
(af + — 2.08) Fy (s, + — (af +08 — 203) Fy (8, p) 
hel 1 (8, P) 3 i+ 3) (8, 
(a8 — a) Fy (8, + 01005 (08 — a4) Fy (8, p) 
mit 4 beliebigen Funktionen F; von s und p. Will man f, und f, in Form einer a 
Reihe nach orthogonalen Funktionen erhalten, so muß man die Funktionen F, 2 
in Potenzreihen entwickeln und die erhaltenen Teilfunktionen orthogonali- 
sieren und normieren. Für die niedrigsten Werte von / ist das von anderer 
Seite längst ausgeführt worden’). Die Durchführung für beliebige / ist uns 
nicht gelungen. Man kann das Verfahren aber so vereinfachen, daß die Ortho- 
gonalisierung bis größenordnungsmäßig ! = 15 ohne viel Mühe gelingt. Auf 
; die Darstellung dieser Methode soll hier aber verzichtet werden. 
A: sd Auf dieselbe Weise kann man die allgemeine Form der Funktionen zu anderen 
ea: Darstellungen finden. Setzt man das Ergebnis in (1’) ein und beachtet (1), 
uinen} so erhält man z. B. für die erste Matrixkomponente mit anderen, von den 
man} obigen verschiedenen Funktionen F, (s, p) für die ea O: 
Ay = F, + of Fy + af Fy +0,00 (038 {Fy + of + af Fe}. (13) 
Ag, und Ag, gehen daraus durch zyklische en hervor. Das läßt 
sich noch vereinfachen. Setzt man in s und p und im ersten Teil von (13) 
3D x? = 1 — 08 — a8, so sieht man, daß F, + af F, + «fF, eine symmetrische 
Funktion S (a3, «%) von a3 und «2 ist, also eine Funktion, für die S (a3, a3) = 
S (x%, «3) gilt. Umgekehrt kann auch jede solche analytische Funktion bei 
ausreichend starker Konvergenz ihrer Potenzreihenentwicklung in die Form 
F, + a? F, + a4 F, gebracht werden. Denn wegen der Symmetrie können die 
Glieder der Potenzreihe so zusammengefaßt und umgeordnet werden, daß nur 
noch Potenzprodukte von a3 a3 und a3 + a3 vorkommen und durch noch- 
malige Umordnung wegen «3 + «3 = 1—a? und a$a3 = s+ af—aj nur 
noch von a? und s. Die Potenzreihenentwicklung hat dann die Form S (a3, «3)= 
5 (a})" - f„(s). Nun gilt weiter «f = p + af — a? s, und damit af = af p + p 
ah: +-at—-a?s—afts, usw. Jede Potenz (x})* mit n > 3 läßt sich somit schritt- 
d= weise in ein Polynom aus s und p verwandeln, in welchem x? nur in 1. und 2. 
t ge: Potenz vorkommt. Setzt man das oben ein, so nimmt S (x3, x3) die verlangte 
Lälfte Form an. 
Natürlich kann man die hier am Beispiel der Gruppe O geschilderte Methode 
tent auch auf die nicht kubischen Kristallgruppen anwenden und auf diese Weise _ 
a die nach der ersten Methode erhaltenen Ergebnisse kontrollieren. Die Ergeb- 
des | nisse, die nach diesen beiden Methoden erhalten werden, sind in den folgenden 
olten Tabellen zusammengestellt. 
Ein 7) F.C. von der Lage. und H.A. Bethe, ar Rev. 71, 612 (1947); D. G. Bell, 
rails | Rev. Mod. Physics 26, 311 (1954). 
E 
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IV. Ergebnisse 


In Tab. 2 sind die nach der zuletzt beschriebenen Methode erhaltenen 
Resultate für alle Kristallgruppen zusammengestellt. Mit großen lateinischen 
Buchstaben bezeichnete Funktionen hängen jeweils ab von den zwei im Kopf 
der Spalte genannten Argumenten. Es ist zu beachten, daß bei manchen 
Kristallgruppen noch Funktionen mit anderen Argumenten oder dieselbe Funk- 


Tabelle 2 
Die allgemeine Form der Tensorkomponenten 4;; für alle Kristallgruppen 


(2453 08)| (03 03) 


C, | | S, 


| (%; Cen [GE %s) 


| 
B+ a, 0, C | B | 


Ay B G, + a, U, G, B(a,, a7) 

Ass C G, + a, U, G, D+oo,F! D 03) 

Ags K G, + a, U, G, K + 0, 0 DL | K K 

Ars M t+ a, U, G, 8 M 

Ags P G, + «, U; U, sP+%0Q %Q P(a3, 03) 
31 a3Q Gy, + a, U, a, U, a, R a, P(a,, 07) 


Cy | (053%) 


Ay B(a, 0%) + 8, 0?) | B+cC i 
Ass B(a3, 03) — 83 03) | B(— ag, — % %&) — Cy C(— ag, — 9%) 
As K+s8,L G, +c U, 
Ais co M + U,+cG, 
Ass 
As, P(— — — Q(— ay, — a, 
D, | (09308) | Di | (af; 03) | Daya 

Au B+ ay | B B(07, 03) + X a3 03) 
Ags D + a4 D B (ag, 05) + & a a3) 
Ag, K + K G (a1 3) 
Aj. x, M + a, a, % N U (a, 
Ass | a, P (ag, 08) + xy Q(a8, a2) 
As R + a, a, T Xy Xy P (xz, 02) + a, a?) 

|» | (a2; a3) | Din | (509) | C,, | (at; a3) 

An | C(af, 02) | Bla}, af) G (044, 03) 
Ag. | Bag, a3) — 0, B(a?, a2 G (— 
Ag | Sı + A, | S, S,+s,A 
An | + a, A, % S, 8, S, + c, S, 
Agg | % P(a§, 03) + (a3, a2) %%Q(08,03) P(ad, a2) + 0,Q (a8, 
Cs | (3 63) | Se | C3) 

Aj, B+3, C+a, D+a, H+c, F+s, H G, +3, U, U, +8 
Ass B+3, C—a, D—a, E—c, F—s, H Gy+s8,U)—«,U,—«,U,—c.G,—38,.G, 
Ags K+s,L | 
Asa +, U, + 8, G, G, 
Ass U,+a,U,+¢,G,+ 8, G, 


A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
A 
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5 4 
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| Cs, | (3 63) | D; | 


(03, Cs) 


Ay B+a,D+ceF 


B+ a38,C + «4, D+ +c, F + H 


Ass B-—a,D—c,F B+ — a, D— a, 0, —c, F — 058, H 
Ags K | K + a, 8, L 
—o,D+s8,F D+ 8, F — H 
Asa P+ 8, T P+ 0,90+¢,R + 0,87 
Ag, x, P—c,T P — Q — 8, T 
Caan | (3; ¢3) |Da |! 


Ay B+s, C+a, D+o, E+c,F+s, H 
Ags B+8, C—a, D—a, E—c, F—s, H 
Ags K+ 

—c,H 
Ass 8, T) 
X3(%, P— + 8, R — T) 


B+oa,D+¢,F 
| 
K 


B+coD+aF 
B-6D-«F 
K 


| 
P — asc, T 


8, D—s8,F 
+ 38,7 


Ce | 


Cg) 


| Con | ( 35%) 


Ay B+8, C+c, D+s8, E+c, F+s, H 
Ago B+8, C—c, D—s, E—c, F—s, H 
Aggy K+a,L 

Ao 8, 

Ags 

Ay, —% P+a,90+s8,R—c,T 


| 


B+s, C+c, D+s, E+c, F+s, H 
B+8, C—c, D—s, E—c, F—s, H 
K+8,L 


P+a,9+c,R+ 7) 
a, R—c; T) 


Da | (035 Cg) 


Dz | 


(03; Cg) 


Ay B+-%3¢3C + «30, 
Ag. | B+a3¢,;C—a,«, D—c,E—c,F—o,¢,H 
K + a3 ¢, L 

— %,%,D+ 8s, —s, F + «38; H 
Ags Xe te P+ 38,90+8,R + 0,8, T 


Ag, P—c, 


Qic, R—a,¢,T 


B+038,C +c, 
B+a38,C—c, D—a,8,E—c,F—o,8,H 
K + a, 8, L 
8, D — a,c, E —8,F + 
% 4,8, 7 
P—a,Q +8, R— age; 


Don | 


(835 Cg) 


Ay, B+6D+«4F 


Ass 
Ass K 
Asa 
7 
Ags P + a3 8, T 
As P—age, T 


T | (035 03) |T'n| (033 03)| O | (085 08) 


(a8; 02) |O,| a8) 


Ay | B+ a, 0,0,C | B 
Ag | a,a,D 


Ago, Az; wie A,, mit zyklischer Vertauschung der Variablen 
A3,, Ajg Wie A, mit zyklischer Vertauschung der Variablen 


8, +9,00 Ay 
%0% 8, + a, A, 


| Te | 


S, + a S, 
| S, + a, S, 
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tion mehrfach, aber mit veränderten Argumenten, vorkommen können. In 


diesen Fällen sind die Argumente direkt angegeben. Die mit fetten Typen 
bezeichneten Funktionen haben noch eine besondere Bedeutung: es ist S 
eine symmetrische, A eine antisymmetrische, @ eine gerade und U eine un- 
gerade Funktion der Argumente: 


S (x, y) = S (y, 2), G (x, y) = G (— x, — y) 


Cn = + 044)" + — 
5 (lo + og)" — (044 — %)”). 


Die Tab. 3 bedarf einer etwas ausführlicheren Erläuterung. 
sultate der anfangs beschriebenen Methode für die 26 uniaxialen Kristall- 
klassen zusammen. Spalte 3 enthält die Bezeichnungen der irreduziblen 
Darstellungen nach Eyring, Walter und Kimball®). In Spalte 4 sind die 
der Kristallsymmetrie angepaßten Linearkombinationen der ß,ß, den irre- 
duziblen Darstellungen zugeordnet, während in den Spalten 5 und 6 die zu den 
Darstellungen gehörenden Kugelflächenfunktionen angegeben sind. Zu be- 
achten ist, daß bei der Bezeichnung der irreduziblen Darstellungen nach 
Eyring, Walter und Kimball der Buchstabe E sowohl für eindimensionale 
wie für zweidimensionale Darstellungen benutzt wird. Dabei sind die zwei- 
dimensionalen Darstellungen nur bis auf Äquivalenzen gegeben. In unserem 
Falle sind jedoch die Darstellungsmatrizen durch die Wahl der Linearkom- 
binationen ß; ß, bestimmt. Die zugehörigen Kugelflächenfunktionen spannen 
genau dieselbe Darstellung auf. 

Ein Beispiel zeigt am schnellsten, wie man verfahren muß. Für die Gruppe 
C, hat die Magnetostriktion die Form: 


A = (Bi +B 3) pi + B3 yd + (67 — 63) pF + B, yt + (8, + Be) Bs ys 
+ (By — By) By 


5 Yim »= vi m 

Y = ye m > Vt m 
m in+2 

m=4n+l 


Da die Magnetostriktion A reell ist, müssen in diesem Fall y, bis y, reell 
sein und y, = ys- Bei geeigneter Umordnung der Glieder kann man also 
y, bis y, sowie y$ + ys und i (y¥ — yg) als Reihe mit reellen Koeffizienten 
und reellen Kugelflächenfunktionen schreiben. 

Bei den Kristallgruppen mit einer zweidimensionalen Darstellung müssen 
die Funktionen zur ersten und zur zweiten Zeile die gleichen Entwicklungs- 


8) H. Eyring, J. Walter undG.E. 
1949. 
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Die in der Tabelle erscheinenden Größen c, und s, sind wie folgt definiert: 
=> 
| 
Ay 
Ad 2x 
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Po Kimball, Quantum chemistry. Wiley, New York 
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koeffizienten haben. Z. B. wird in der Kristallgruppe D, die zweidimensionale 
Darstellung E durch (8, 83) aufgespannt. Das reelle Funk. 
tionenpaar ist (y;; y,) mit 


l=2k 1=2k+1 


1=2k+1 
9 
2n+1 m=2n+1 
Ae 


V. Anwendungsbeispiele 


_ Als erstes Anwendungsbeispiel betrachten wir die Magnetostriktion eines 
hexagonalen Kristalls der Symmetrie D,,. Setzt man die in Tab. 2 enthaltenen 
Ausdrücke in (1) ein, so erhält man nach geeigneter Zusammenfassung den 
folgenden Ausdruck mit 6 unbestimmten Funktionen von a} und c, = aj 


— 15 af a3 + 15 af af — af 


A B;) = (Bi + B+ {BT — PB) — x8) + 4B, By D 

{(Bi — (at + a8 — 60703) — 8 B, (af — F + PRK 

+ 2 (By 04 + By Bs P + 2 (By — Be 85) T. 

Bei der Kerala dieser Formel ist zu beachten, daß man üblicher- 
weise die magnetostriktive Dehnung von demjenigen Zustand aus mißt, 
in welchem die Magnetisierung gleich mäßig auf die leichten Richtungen ver- 
teilt ist. Im Kobalt ist dann 9 = 1, a, =a, = (0. Wählt man auch hier 
diesen Nullpunkt, so ist in B und K die Konstante so zu wählen, daß diese 
Funktionen für x = 1 und c,= 0 verschwinden. Entwickelt man die 6 
Funktionen B, D, F, K, P, T nach Potenzen ihrer Argumente und berück- 
sichtigt überall die Glieder bis zu den 4. Potenzen von «,, so erhält man die 
von Mason bereits abgeleitete Formel mit 9 Konstanten in anderer Zusammen- 
fassung wieder. Die obige Formel erlaubt es, fast ohne Rechenarbeit die Ent- 
wicklung beliebig fortzusetzen, was bei der Rechenmethode von Mason 
bald zu sehr unübersichtlichen Berechnungen führen wiirde. Berücksichtigt 
man überall noch die 6. Potenzen der «,, so ist zu setzen 


B=b(1-o5) +b,(1— 03% +b,(1— 05)’ + by cg, 

K=k, (1—§) + & (1 — 08)? + (1 — + cy, 
+ 403 + dy af, F=f,+ f, «3, Per 

T=t. 

Dann treten also 17 Entwicklungskoeffizienten auf. Erst bei Berück- 
sichtigung dieser Glieder hängt die Magnetostriktion nicht mehr allein von den 
Winkeln zwischen der Meßrichtung, der Magnetisierungsrichtung und der 
hexagonalen Achse ab, sondern außerdem von den Winkeln zwischen den 
zuerst genannten Richtungen und den zweizähligen Achsen senkrecht zur 
hexagonalen Achse. In der gleichen Weise wie in diesem Beispiel liefern die 


Tabellen allgemeine Aussagen über die Abhängigkeit der Magnetostriktion, 
der und anderer tensorieller von der Ma- 
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achten, daß hier noch eine Symmetrieforderung hinzukommt, die nicht aus 
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Die Anwendbarkeit der Ergebnisse ist aber damit nicht erschöpft; sie 
gelten ja für jeden Tensor, dessen Komponenten Funktionen einer Richtung 
sind. Ein weiteres Beispiel für einen solchen Tensor liefert die Theorie der 
kleinen Schwingungen eines Gitters. In ihr tritt ein Tensor ® auf, welcher den 
Zusammenhang zwischen der Kraft % auf ein erstes Gitteratom und der Ver- 
rückung irgendeines zweiten um den Vektor $ aus seiner Ruhelage bei unver- 
schobenen anderen Atomen angibt?): 


— 


Für ein bestimmtes erstes Atom ist dieser Tensor ® noch eine Funktion 
des Ortsvektors r von dem ersten zu dem zweiten Atom. Er ist invariant 
gegen alle Symmetrieoperationen des Gitters, die das erste Atom in seiner 
Lage lassen. Zum Unterschied von den früheren Betrachtungen ist hier ® 
nicht nur von der Richtung des Vektors r, sondern auch noch von seinem Be- 
trage abhängig. Wenn man aber nur die Tensoren für diejenigen anderen Atome 
vergleicht, welche durch die genannten Symmetrieoperationen auseinander 
hervorgehen, spielt das keine Rolle, weil |r| dabei unverändert bleibt. Im 
einfach-kubischen, kubisch-flächenzentrierten und kubisch-raumzentrierten 
Gitter hat die Umgebung jedes Atoms die Symmetrie O,. Die Tensoren für 
alle zweiten Atome im gleichen Abstand von dem ersten, haben also die Form 


Sy (03, 09) Sy (HF, 03) Sy (aF, 03) 
Sy (af, S, (a3, af) Sy (x3, 3) 
X] S, (af, a3) Ag S, (x3, a) S, (a7, 3) 


Betrachtet man ein einfach-kubisches Gitter und die nächsten Nachbarn, 
so gibt es ein Paar von Nachbarn mit «, = + 1, % = & = 0 und zwei weitere 
Paare, deren Richtungskosinus durch zyklische Vertauschung erhalten werden. 
Schreibt man S, (0,0) = a@ und S, (1, 0) = S, (0, 1) = b, so lautet die Matrix _ 
fir diese 3 Nachbarnpaare 5 


ado 5600 600 
1.Paar [05 3. Paar{[ 0060 
005 006 00a 


Die bei Leibfried angegebenen Tabellen dieser Matrizen kann man auf 
diese Weise praktisch ohne Rechnung hinschreiben. Es ist nur noch zu be- 


der hier betrachteten Symmetrie des Gitters folgt. Da die Komponenten von 
® in diesem Fall zweite Ableitungen einer potentiellen Energie sind, muß für 
Atompaare, die durch eine Translation des Gitters ineinander überführt 
werden können, ein Vorzeichenwechsel der a, einer Spiegelung der Matrix an 
der Hauptdiagonale gleichwertig sein. 


°) G. Leibfried, in Hdb. der Physik, herausgeg. von S. Flügge, Bd. VII, 1, S. 104ff. 


Gießen, Institut für theoretische Physik der Justus-Liebig-Universität. 


Bei der Redaktion eingegangen am 25. Juni 1959. 
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"Zur Theorie der Wechselschichten 
aus schwachabsorbierenden Substanzen und ihre 
Verwendung als Interferometerspiegel 


Fon Koppelmann) 


Es wird mit Hilfe der Matrizenmethode die Absorption von A/4- und 4/2. 
Wechselschichten aus schwachabsorbierenden, dielektrischen Substanzen 
berechnet. Ferner wird der Einfluß der Absorption bei Verwendung dieser 


Schichten als Interferometerspiegel untersucht. 


> 
= 


1. Einführung der Matrizenmethode 


Auf ein System von p planparallelen, homogenen, isotropen Schichten 
falle eine monochromatische, ebene elektromagnetische Welle senkrecht 
auf (s. Abb. 1). Innerhalb jeder Schicht ergibt sich dann eine einfallende und 
eine reflektierte Welle mit der Amplitude € bzw. kt. Mit Hilfe der Maxwell- 
Gleichungen und unter Benutzung der Stetigkeitsbedingungen an den Grenz- 
flächen ergibt sich ein lineares Gleichungssystem für die (komplexen) Am- 


plituden € bzw. i des elektrischen Vektors der 
| Welle innerhalb der einzelnen Schichten?). 
(0) 
et ie — on, Im, — 
Hierbei ist n, der komplexe Brechungsindex des 
(p-1) 
Mediums s 
(p) =n, —ik, (2) 
(ptt) und 
Abb. 1. System von p Teil- 7 ee 2a n, d,. (3 
schichten 


I) Aus der Dissertation des ag ee 1959. 

2) Eine ausführliche Darstellung und Ableitung s. H. Mayer, Physik dünner Schichten 
(Teil 1), Stuttgart 1950. Die (1) entsprechende Formel (8.171) bezieht sich dort jedoch 
auf die magnetische Komponente des elektromagnetischen Feldes und weicht daher in 
geringem Maße ab (Faktor n,/n,+,)- 


6. Koppe 
wobei d 
ist die | 
reflek tie 
der Gre 
Fühı 
ein 

ts 
folgt 


Die- 
zehung 
und nac 


Durch : 


wobei 


und dar 


Also sin 


Die : 
keit sin 
j 


p 


/ 


» 


Wir 


system 


von (di 


4) p ‘ 
tteffende 


| 
= 


6.Koppelmann: Theorie der Wechselschichten aus schwachabsorbierenden Substanzen 389 


wobei d, die Schichtdicke und A die Wellenlänge (im Vakuum) ist. €, bzw. R 
it die komplexe Amplitude des elektrischen Vektors der einfallenden bzw. 
reflektierten Welle im Medium (s) und zwar gemessen in unmittelbarer Nähe 
der Grenzfläche (s — 1); (s). 

Führt man nun die Fresnelkoeffizienten für die einzelnen Grenzflächen 
an 
— 2n, 


Tsti = 
so folgt 


wg 


= + ei Resa. 
stl st+1 
Die: Gl. (5) stellen einen Satz von Rekursionsformeln dar, die eine Be- 
tichung zwischen den Amplituden Rig, E, und €, herstellen (6) = 0; R,+] = 0) 
A/2-Jund nach Abelés*) auch in Matrizenform geschrieben werden können: 


nzen id 
Durch sukzessives Einsetzen unter Beachtung von R,+; = 0 erhält man 
(7) 
und 
‚ejj. [und damit nach (7) 
1 
&=a:6,4ı Ryo = (9) 
Am- ty to+1 t, toss 
; derf4lso sind die gesuchten 


Die zugehörigen Intensitätsverhältnisse der Reflexion und der Durchlässig- 
keit sind dann®) 


| 


fe 0 
* 
T — Np+1 E,+1 1 t, t*. -% 2 
| & My aa*s 
des 2 
2) 2. A/4-Wechselschichtsystem aus . 7 
( schwachabsorbierenden Substanzen 
(3) Wir betrachten jetzt ein spezielles Schichten- % 
‘' [system (Abb. 2), das aus einer ungeraden Anzahl = 
von (dielektrischen) Schichten mit abwechselnd “ts 7 


rn *) F. Abelös, Ann. Phys., Paris, 8, 504 (1948). 
edoch 4) p als Index links unten deutet an, daß die be- Abb. 2. Wechselschichten- 
1er teffende Größe für eine p-fache Schicht gilt. system 


| 
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hohem und niedrigem Brechungsindex (n,; 115) besteht. Die oberste Schicht 
habe einen hohen Brechungsindex n,; die Brechzahl des vorderen Mediums 
sei %,, die des Schichtträgers n,. — Für derartige Systeme mit absorptions. 
freien Schichten (also n, und n, reell) und den optischen Schichtdicken 
n, dy = Nd, = A/4 lassen sich die Matrizenprodukte (8) relativ einfach aus. 
rechnen, da alle Phasen 6, = 2/2 werden (außer ö, = 0), wie bereits Heavens 
gezeigt hat. 


Nehmen wir nun an, daß zwar ebenfalls ein A/4-Wechselschichtensystem 
vorliegen möge, in dem aber die dielektrischen Schichten eine gewisse Absorp- 
tion besitzen mögen, so sind die Größen n, und damit auch r, und 6, komplex 
und damit die explizite Berechnung des Reflexionsvermögens und der Durch- 
lässigkeit auch nach der Matrizenmethode sehr kompliziert. Im folgenden 
soll nun aber gezeigt werden, daß man zu relativ einfachen Näherungsformeln 
kommt, wenn man ansetzt, daß die Absorption in den Schichten nur klein 
sein möge (k, < n,), was für wickitanteiiische Interferenzschichten im all- 


gemeinen zutrifft. 
Die Phase 6, ist nach (2; 3) SV 


n,d,— k, d, fir s=1;2;...p und 6=0. (12) 
Für 4/4-Schichten ergibt sich mit der Abkürzung e, oe Rr 
(13) 


3 2a 


. 
Fi ür schwachabsorbierende Substanzen (k, < n,) st, => 


-— <1: daher 
n, 
kann man näherungsweise ansetzen 
; 


? 
ußerdem nac ch (4) = —t3 = =—tr,; damit folgt aus (8) 


+ &) (U + & 0); —(1+4)9 y= ( te + 


wobei gesetzt ist 
u=1+% w=1-1 v= wiki (15) 


Durch Ausrechnen der Matrizenprodukte, wobei immer die in ¢ quadra- 
tischen und höheren Glieder gegen die linearen vernachlässigt werden, erhält 
man (nach einer hier nicht wiedergegebenen Zwischenrechnung) die folgenden 


0.8. Heavens, J. opt. Soc. Amer. 44, 371 (1954) 
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shicht} Formeln (17), wobei außer den obigen Abkürzungen benutzt wird 


1 151 tr = A 1 tr = € (16) 

— -.= = ‘ 
at =A +46 


stem p= 


\SOrp- | c(i)? = +e,wD (17) 
nplex _ > = Au? + 2Buv + Av? + [3u? mE [2u] 
= — Bu? — 2A uv — + D [3u? + v7] + D [2u] 
"meln 

klein + 


‘? Für das allgemeine Glied dieser Reihe ergibt sich nach Durchrechnung mit 


(15) (16): 


| (Mp + My) p+1 1 
Girne (18) 

1 (tg — My) ine 

+ ny — Ng 

laher wobei 
1 
(My + Ty) (1, + Ny) an, 2n, 


Aus den Größen a und c ergibt sich nun nach Formel (11) das Reflexions- 
und Durchlässigkeitsvermögen ,R und „T des gesamten Schichtensystems. 
sw Die Brechzahlen n = n — ik in der Formel (18) sind im allgemeinen kom- 
plex, wobei jedoch k <1 ist. Die Größen a(i)-? und c(i)-? sind also ebenfalls 
komplex, wobei der Imaginärteil nur erste und höhere Potenzen von k ent- 
hält. Bei der Bildung von a a* und c c* trägt der Imaginärteil also nur mit 
quadratischen und höheren Gliedern von k bei, die nach Voraussetzung zu 


tz! vernachlässigen sind. Es genügt demnach, in die Formeln (18) nur die realen : 
Anteile der Brechzahlen einzusetzen. 3 
) 
). 3. Die Verteilung der Absorption auf Durchlässigkeit und Reflexion a 
Infolge der Absorption wird einerseits die Reflexion „R und andererseits 
die Durchlässigkeit „7 um einen bestimmten Betrag A,R bzw. 4,T gegen- a 
(19) über der eines entsprechenden, absorptionsfreien Schichtensystems „R, bzw. ; 
pl) zurückgehen 
hält ot = ,R, —A,R =,17,—A,T. (20) 
den Man kann diese Tatsache so auffassen, als ob sich die Gesamtabsorption „A 


auf Reflexion und Durchlässigkeit „verteilt“. 
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„R und „T berechnen sich nach Formel (11), wobei sich entsprechend (18) p> 00°) 
für a und c (da p ungerade, praktisch rein imaginäre) Größen von der Form fjweils e 
ıngenäh 
3. Di 
ergeben, wobei a,;c, Konstante und a,, c, linear in e, und ¢, und damit in  |Shichte 
k, und k, sind; also ist immer a,; c, <<a; Cy. a, und c, sind dabei die ent-  |indert s 


a=a,+a, C=C) 


sprechenden Werte im absorptionsfreien Fall. Diese 
Für die Durchlässigkeit ergibt sich daher®) sondern 
den For 
1 —. n + 2 P . 
= (a, + a,)* No = + (1 + N „Io ( | 
0 € == I, 0 > 0 
und für die Reflexion Bei Betı 
e,\? 1 + 2¢,/c 2c, 2a ich auf 
im a, +a, (1 + 2a Ja) Co dy 
also 
2a 2a 2c ü 
A,T=,%— A,R=,R, > >). (24) 


Die hier abgeleiteten Formeln sollen nun an einem Beispiel erläutert werden. 
Wir betrachten speziell ein hochreflektierendes A/4-Schichtensystem mit 
N, = 2,30; ng = 1,30 und np = 1,50 — entsprechend einem System von 
ZnS-Kryolith-Schichten, das sich auf einem Glasträger befindet. Der Ein- also 
fachheit halber wollen wir die niedrigbrechende Substanz als absorptionsfrei 
annehmen. Die Ergebnisse sind in der Tab. 1 zusammengestellt. 


wobei 


Analog 
Tabelle 1 b\ 
b/ 
m A i I A R 
p ols „Ro 4,1 A,R | „A | 
1 0,688 0,312 | 0,798 24 0,230-k | 1,028-k | 1,158-k | 0,738-k in 
3 0,305 0,695 | 0,548 -k| 0,949: k 1,497. k | 1,798. k 1,365 - k 
5 0,109 0,891 0,223 -k| 1,439-k | 1,662-k | 2,050-k | 1,616-& man scl 
7 0,0363 0,9637 0,0775 k| 1,647. k 1,724-k | 2,1388-k | 1,708-k 
9 0,0118 0,9882 0,0256 -k| 1,712-k |1,738-k | 2,165-k | 1,732-k a (— 
1 0,0037 0,9963 | 0,0081 -k| 1,734-k |'1,742-k | 2,.179-k| 1,741-k 
00 0 | 41 0 | 1,744-k | 1,744 1,744 k 
Aus der Tabelle entnehmen wir: Nach d 
1. Bei den hochreflektierenden Mehrfachschichten (etwa ab p = 5 Teil- er 
schichten) geht die Absorption praktisch vollständig auf Kosten der Reflexion. schichte 
2. Die Gesamtabsorption „A wächst im Bereich hoher Reflexionen nur ) D 
noch wenig mit steigender Schichtenzahl und strebt einem Grenzwert zu für 14-Weel 
fallende 
*) In den an sich komplexen Fresnel-Koeffizienten t, kann man bei der Berechnung Abeorpti 
von „T hier, aus gleichen Gründen wie oben, den Imaginärteil von n = n — ik vernach- » ander 
lässigen und mit den reellen Größen r, = rechnen. 7 Di 


> 
e 
4 
7 
° 


§,Koppelmann: Theorie der Wechselschichten aus schwachabsorbierenden Substanzen 393 


oo”). Auch die Verhältnisse A,7/,T, und A,R/,R, nähern sich für p > 
weils einem Grenzwert. Bereits für p= 5 werden diese Grenzwerte gut 
ıngenähert erreicht. 

3. Die Gesamtabsorption „4 eines Schichtensystems mit unendlich vielen 
Schichten ist nur knapp doppelt so groß wie die einer Einfachschicht; „A 
indert sich also überhaupt nicht sehr stark mit p. 

Diese Aussagen sind nicht auf das spezielle Zahlenbeispiel beschränkt, 
wndern gelten allgemein für hochreflektierende A/4-Wechselschichten. Aus 
den Formeln (18), (24) ergeben sich für 2— oo (wenn man beachtet, daß 
n/n)” —> 0 geht) die sehr einfachen Formeln für die Grenzwerte®). 


Aw? _, (hy + k,) + k,). (25) 
cof, 
Bei Betrachtung hochreflektierender ON A = 5) kann man 
ich auf diese Formel als gute Näherung beschränken. 
4. Die Absorption eines A/2-Wechselschichtensstems — 


Für A/2-Schichten haben wir in die Formel (8) einzusetzen 


2 
— d,, (36) 
N, d, > = 7 k, d, = Ns s’ oF 
Analog zu Formel (14) ergibt sich dann die Matrizengleichung (27) AR 
cy rt, 1 (1-e) (1-6) 
-( (1 + &) 3 + (27) 


Rechnet man nun in ganz analoger Weise die Matrizenprodukte aus, so erhält 
man schließlich für p Schichten mit den Abkürzungen (15) (16) ~ i 


(28) 


Nach den Formeln (21); (24) kann nun leicht „A; A,R und A,T berechnet 
werden. 

Als praktisches Beispiel wollen wir wieder das ZnS-Kryolith-Wechsel- 
schichtensystem betrachten, in dem nur die ZnS-Schichten absorbieren sollen. 


?) Das Auftreten dieses Grenzwertes läßt sich auch anschaulich versteben: Von einem 
i4-Wechselschichtensystem mit zahlreichen Einzelschichten wird der Hauptanteil des ein- 
illenden Lichtes bereits in den obersten Teilschichten reflektiert; dort findet auch die 
Absorption im wesentlichen statt. Fügt man nun noch ein weiteres Schichtenpaar hinzu, 

% ändert sich der Reflexionsvorgang und damit die Absorption nur geringfügig. 

8) Diese Formeln stimmen mit den von P. Giacomo, Rev. Opt. (theor. instrum.) 

%, 317 (1956) auf anderem Wege abgeleiteten Formeln überein. ee er? or 
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Ri ie Die Ergebnisse sind in der Tab. 2 zusammengestellt. — Eigenartigerweig 
er eben sich jetzt für 4,R negative Werte, d.h. das Reflexionsvermöge 
8 ge 
steigt infolge der Absorption etwas an?). RN 
ait 


1 3,09-k | —0,299-k 286-k 

| | | 618-k  —0,459-k 5,72. k 

5 0,96 0,04 9,27-k | —0,689-k 8,58 -k 
7 | | 1236:k | —0919-k | 

4 15,45-k | —1,160-k | 14,80-k 

, kar Vergleicht man nun die Absorption der A/2-Mehrfachschichten mit der der 


iro At Mehrfachschichten (s. Tab. 1), so erkennt man, daß die Absorption bei 

gleicher Schichtenzahl p im A/4-Fall ganz wesentlich geringer ist als im Aß. 
m Fall. Insbesondere steigt sie im A/2- Fall mit wachsendem p sehr stark an (nach 

7 _ einer arithmetischen Reihe), während im //4-Fall schon sehr schnell ein relativ 


niedriger Grenzwert praktisch erreicht wird. 


x; Man erkennt schon aus der Tatsache der verschiedenen Absorptionen von 
‘ta 4/4- und A/2-Schichtensystemen, daß man bei einer dielektrischen Mehrfach. 
_ sehicht mit schwachabsorbierenden Teilschichten nicht mit einer von d 
Wellenlänge unabhängigen Absorption rechnen darf!P); vielmehr hängt die 
_ Absorption stark davon ab, ob die Mehrfachschicht für die betreffende Wellen- 
länge gerade hochreflektierend ist oder stark durchlässig. 


Die Berechnung der Absorption einer Wechselschicht für eine beliebige 
Wellenlänge A führt zu sehr komplizierten Formeln. Wir haben hier die Absorp- 
= tion für die beiden Grenzfälle der A/4- und A/2-Wechselschicht berechnet. inf 
allgemeinen Falle dürfte die Absorption zwischen diesen Grehzfällen liegen. 
ni Man kann übrigens leicht zeigen, daß kleine Abweichungen von der A/4- bzw. 4/2 
Bedingung keinen Einfluß auf das Absorptionsvermögen haben: In diesem Falle ist nam, 
lich statt (13) bzw. (26) einzusetzen 6, = 2/2 + y—ie, bzw. , =n+y-—ie, Wit 
erhalten also die gleichen Ergebnisse wie im A/4- bzw. 1)2- Fall, wobei nur das reelle ¢ 
zu ersetzen ist durch das komplexe (¢,— iy). Da nun e nur linear inaund ca 


De Pfs und nur a a* und c c* in die weitere Rechnung eingeht, tritt y nur als y? auf und ist gegen- 
Sr über den in e linearen Gliedern zu vernachlässigen. 
| 
os 5. Verwendung schwachabsorbierender A/4-Wechselschichten 

als Interferometerspiegel 
= Für ein Pérot-Fabry-Interferometer mit zwei Spiegeln, die jeder f 
AR sich ein Reflexionsvermögen R, eine Durchlässigkeit T und ein Absorptions- 


vermögen A besitzen, ergibt sich bekanntlich eine maximale Durchlässigkeit J, 


nena 9) Betrachtet man zur anschaulichen Deutung dieses Effektes z. B. die zwei an de 
_ Vorder- und Rückseite einer freitragenden A/2-Schicht reflektierten Strahlen 1 und} 
so löschen sich diese im absorptionsfreien Fall vollkommen aus, da sie gleiche Amplituds 
er und entgegengesetzte Phase besitzen. Wenn dagegen die Schicht schwach absorbiert 


R so ist die Amplitude des Strahles 2 geschwächt und es wird bei der Interferenz eine Rett 
Er reflexion übrig bleiben. 
S. Penselin u. A. Steudel, Z. Phys. 142, 21 (1955), 
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innerhalb der Interferenzstreifen und eine Apparatebreite h (d. h. Halbwerts- 
breite bezogen auf den Streifenabstand) entsprechend den folgenden Formeln 


(wobei die Näherung für R ~ 1 gilt). 

Es mögen die beiden Interferometerspiegel nun aus zwei dielektrischen 


4/4-Wechselschichten bestehen, die wir zunächst als absorptionsfrei annehmen 
wollen. Dann ergibt sich für das Reflexionsvermögen aus (11) und (18) 


(30) 


Für höhere Schichtenzahlen ist (n,/n,)? < 1 und man kann dann näherungs- 
weise ansetzen 
jl = = [1 - N, (=) why 1—,R, (>) (31) 
nz 


Ny an \n 


Wie schon Giacomo?) gezeigt hat, wächst also die Schärfe der Interferenz- 
streifen etwa um den Faktor (n,/ns), wenn man in den absorptionsfreien 
Interferometerspiegeln die Schichtenzahl um 2 erhöht. 

Besitzt nun in den A/4-Wechselschichten jede Schicht eine gewisse, kleine 
Absorption, dann nimmt, wie wir gesehen haben, sowohl das Reflexions- 
vermögen R, um einen bestimmten Betrag AR als auch die Durchlässigkeit 7, 
um den Betrag AT ab. Die Absorption ist dann A = AR + AT. 

Betrachten wir schlieBlich die Frage, wie sich Apparatebreite und Durch- 
lässigkeit eines PFI verhalten, wenn man die Schichtenzahl in den Spiegel- 
systemen erhöht. Geht sie nach oo, so nähert sich sowohl A als auch AR 
nach (25) einem gemeinsamen Grenzwert 4,.R = „A, während die Durch- 
lässigkeit 7’ gegen Null geht. Daraus folgt, daß mit steigender Schichten- 


zahl /,>0 und h> od geht. Es läßt sich also die Halbwertsbreite 


auch bei Verwendung von Spiegeln aus dielektrischen Mehrfachschichten 
nicht beliebig herabsetzen, sondern nur bis zu einem gewissen Grenzwert 


ak = — oA. Dieser beträgt nach (25) 
Zu 


Betrachtet man eine endliche Schichtenzahl, so erweist sich die Größe 
z= (33) 


— also das Verhältnis der Durchlässigkeit der p-fachen, absorptionsfreien 


Schicht zum Absorptionsvermögen im Grenzfallp > oo — vonentscheidender 


Bedeutung für die Abweichungen der Interferometereigenschaften (Halb- 


wertsbreite „A und maximale Durchlässigkeit „/,„) von denen einesentsprechen- 


den, absorptionsfreien Interferometers („Ag; „/mo = 1). Berücksichtigt man, 


daß bei hochreflektierenden Systemen (schon für p = 5) näherungsweise gilt _ 


und 
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so folgt als Näherung nach (29) Räte 
2 
In Abb. 3 sind diese Größen als Funktion von x aufgezeichnet. Solange x Me 
wesentlich größer als 1 ist, weichen die Eigenschaften des realen PF-Inter- 
ni ferometers von denen eines absorp- 
Pr tionsfreien nur unwesentlich ab. 
Fr % Geht jedoch x gegen 1, oder wird 
ah = Le es sogar noch kleiner, — d. h. liegt 
—— pln 60 die Durchlässigkeit des Spiegel- 
wi schichtensystems in der gleichen 
- + 40 Größenordnung wie seine Absorp- ¥ 
Va tion —, so nimmt die Durchlassig- deren 
2r- 20 keit sehr schnell ab und die Halb- bestir 
0 0 wird 
5 persed Herrn Prof. Dr. K. Krebs bin Verd 
Abb. 3. Die Anderung der Halbwertsbreite h ich fiir sein stetes förderndes Inter- auf A 
und der maximalen Durchlässigkeit Im für esse an dieser Arbeit und für zahl- : 
ein Pérot-Fabry- Interferometer als Funk- unte) 
tion der Größe x = ,Ty/ooA reiche Diskussionen sehr zu Dank 
verpflichtet. Der Deutschen For- 
schungsgemeinschaft danke ich für die Gewährung einer Sachbeihilfe, die mir | 
die Durchführung dieser Arbeit ermöglichte. we 
dem 
Berlin-Charlottenburg, II. Physikalisches Institut der Technischen 
Universitat Berlin. 
und 
Bei der Redaktion eingegangen am 15. Juni 1959. die § 
me [3 


) 
es. Messung kleiner Schichtdickenunterschiede in dünnen 
rp- 
Inhaltsübersicht 
rp- Es wird die Schichtdickenverteilung in dünnen, dielektrischen Schichten, 
sig- deren Dicke sich von Ort zu Ort nur wenig ändert, mit hoher Genauigkeit 


Ib- bestimmt, indem das Reflexionsvermégen in Abhangigkeit vom Schichtort 
gemessen und auf Dickenänderungen umgerechnet wird. Die MeBanordnung 
wird beschrieben, das Ergebnis der Messungen für eine Kryolith-Aufdampf- 
schicht angegeben und daraus die Dampfstrahlcharakteristik der benutzten 


bin Verdampfungsöfen bestimmt. Ferner wird die Anwendung des Verfahrens 
er auf Mehrfachschichten sowie der Einfluß der Dispersion der Schichtsubstanzen 
ant untersucht. 

ink 

or- 1. Die Grundlagen des Verfahrens 

mir 


Das Reflexionsvermégen R einer absorptionsfreien, einfachen, homogenen _ 
Interferenzschicht hängt bekanntlich bei senkrechtem Lichteinfall nur von 
dem Argument 
(1) 
und der Brechzahl der Schicht » und des Schichttragers n, ab; (d ist dabei 
die Schichtdicke und 4 = + die Wellenlänge). In der Umgebung der Stellen, 


an denen die Bedingung 


3 5 
bzw. nd=— (2) 


dem in einem ziemlich großen Bereich in guter Näherung linear (s. Abb. 1). 
Ne ley j lese j ilt ; 1 
Für einen Wert ®, in diesem Bereich gilt also!) don 


dR 
R (®, + A®) = Ry + 7g 4: (3) 
wobei = praktisch eine Konstante ist. 
Wird nun bei konstanter Wellenzahl », die Schichtdicke d um Ad verändert, 
so gilt 
dR 


IR= 19419 = -4anv, Ad, 


1) Die gleichen Überlegungen wie für die Reflexion R(®) gelten natürlich auch für 


die Durchlässigkeit 7 (®). Im Falle reflexionsvermindernder Schichten (s. u.) ist es aber E 
meßtechnisch wesentlich günstiger, Reflexionsmessungen zu benutzen. En 


erfüllt ist, ist die Änderung von R als Funktion von © relativ stark und auBer- 


Annalen der Physik. 7. Folge. Band 5. 1960 ee MER 


fe 


Cz mit C= dy = const. 
Wird andererseits bei konstanter Schichtdicke ‘do die Wellenzahl » um Av ge- 
ändert, so gilt entsprechend 
Ry % 
wobei der Einfluß der ee n (v) zunächst vernachlässigt sei. 

Um nun nach Gl. (4) die Schichtdickenverteilung als Funktion des Ortes 
zu bestimmen, wird mit Hilfe eines schmalen Lichtbündels monochromatischen 
R Lichtes (v9) die Änderung des Re- 
 flexionsvermögens AR als Funktion 

es Ortes (x; y) gemessen gegen- 
über dem Reflexionsvermögen® R, 


j 


in einem Punkt x = 0;¥=0 mit 


é : der Schichtdicke d,. Man wählt 
ry 
0 dabei die Wellenlänge A, = so, 
0 T % en In 9 
’ daB man sich in einem der oben 


Abb. 1. Reflexionsvermégen R einer refle- 


xionsvermindernden Schicht (n = 1,3) als ys nnten linearen Bereiche be 


de findet. Nach (4) besteht dann Pro- 
Funktion von ® = znd portionalität zwischen AR (x; y) und 
Ad(x; y). Die Konstante C, die zur 
Auswertung bekannt sein muß, wird nach GI. (5) durch eine zweite Messung 
bestimmt, bei der man an der Stelle «= 0; y=0 (d.h. d=d,) die Re- 

flexion R als Funktion der Wellenzahl » mißt. 

Um schließlich aus der relativen Dickenänderung Ad/d, den absoluten 
Wert Ad zu erhalten, kann man. d, aus den Extremwerten der Kurve R (») 
ermitteln: Für Lage »,, bzw. v,, und Größe R, bzw. R,, der Extremwerte 
gilt bekanntlich 


2k+1 1 2k pire 
n dy = n dy = ; 2, 


4 Ym 4 vy vr (6) 
n? — n, N„\2 N, — Np \2 
Re = (57 r) bzw. Ry, =(— 
n® + ny Np N; + Nr 
Mit dem Brechungsindex n, des Trägers und n, = 1 (Luft) läßt sich hieraus 
leicht (n dy), n und damit d, ausrechnen. 


2. Durchführung der Messung 
Die von der Schicht reflektierte Lichtintensität wurde in folgender An- 
ordnung gemessen: Der in der Höhe stark begrenzte Monochromator-Aus- 
trittspalt wurde mit einem Hohlspiegel auf die zu untersuchende Schicht ab- 
gebildet (Größe des ausgeleuchteten Bereiches 2 x 1 mm?). Das unter einem 
spitzen Winkel reflektierte Licht wurde dann mit einem zweiten Spiegel auf 
einem Multiplier gesammelt, vor dem sich eine Mattscheibe befand?). Die 
Parallelverschiebung des Schichtträgers um meßbare Beträge x und y erfolgte 


2) Diese Anordnung ist relativ unempfindlich gegen kleine Dejustierungen des zu 
untersuchenden Spiegels während der Messung. 


nit Hilf 
quelle w 
Photostr 

Da es 
des Refl 
sine hohe 
Bs genü; 
I(x; y) 
vergleich 


Im HE 
allein, se 
trägerrü. 


U 


Man hat 


ind ents 


| 


bzw. 


mit 


Die Kon 
Bestimn 
fahren di 
wendet 1 

Die 
mentell 
Reflexio: 
des Träg 
gilt 


Aus den 
R, = 
ty: V 

Die f 
Umrech1 


3) 


~ 
als 
| u 
| — 
7 
|| 
m 
% 
u 


6) 


1S 


G. Koppelmann: Schichtdickenunterschiede nichtmetallischer Schichten 399 


nit Hilfe eines Mikroskop-Kreuztisches. Als relativ gut konstante Licht- 
welle wurde eine Batterie-gespeiste Glühlampe benutzt. Die Anzeige des 
Photostromes erfolgte mit einem Multiflex-Galvanometer. 
Da es hier im wesentlichen nur auf eine genaue Messung der Änderungen 
des Reflexionsvermögens und nicht auf die absoluten Beträge ankommt, ist 
ine hohe Genauigkeit (etwa + 0,1 bis 0,2%, von R,) relativ leicht zu erreichen. 
fs genügt, die an den verschiedenen Orten x; y reflektierte Lichtintensität 
I(x; y) zu messen und mit der Intensität J, an der Stelle «= 0; u = 0 zu 
vergleichen. Es gilt dann 
AR 
9 


Im Experiment wird man allerdings nicht die Reflexion der Schicht R 
illein, sondern eine Überlagerung (o) von R mit der Reflexion der Glas- 
kseite messen?): 


R+! E — R)’ Ry ER er 

Man hat also statt (7) genauer zu setzen > 

I-I__4e 


und e ntspreche ‚hend (4 und (5) 


Die Konstante C* ist auch hier bei beiden MeBreihen die gleiche. Mit ihrer 
Bestimmung nach (5) wird der Beitrag von R, mit eingeeicht, so daB das Ver- 
fahren der Dickenbestimmung nach (4) ohne Umrechnung von o auf R ange- 
wendet werden kann. 

Die spektrale Abhängigkeit des Reflexionsvermögens läßt sich experi- 
mentell recht genau bestimmen, indem das Verhältnis V(v) = o(v)/o, der 
Reflexionen am Ort x = 0; y=0 o(v) und an einer unbeschichteten Stelle 
des Trägers 0, unmittelbar hintereinander gemessen werden. Entsprechend (8) 

Or = Ry 1 =. (11) 

Aus dem — als bekannt vorausgesetzten — Brechungsindex n, des Trägers 
I" Ry = (np — 1)?/(np + 1)?, also auch 0, und damit schließlich o (vy) = 
tp: V (v) berechenbar. 

Die fiir die Bestimmung von d, erforderlichen Größen R,, werden durch 
Umrechnung der Extremwerte o,, der Kurve o(v) nach 


ee 


G. Koppelmann u. K. Krebs, Z. Physik, im Druck. 1% ag 
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= 3. Ergebnisse 


Als ein Beispiel sollen die Messungen an einer Kryolithschicht wieder. 
gegeben werden, die auf eine Glasplatte von 6 cm Durchmesser in einer Ent. 
fernung 22 cm vom Ofen im Vakuum aufgedampft worden war. Es wird dabei 
die Schicht in der Mitte etwas 
dicker sein als am Rande. 
Zunächst wurde die spektrale 
Verteilung des Reflexionsver- 
mögens o (vy) in der Mitte der 
Schicht gemessen. Aus Lage und 
Größe der Extremwerte folgt ent- 
sprechend (6); (11) eine A/4-Dicke 
für v = 1,69. 10% em-!, also nd, = 
1480 Ä, und eine Brechzahl von 
etwa n = 1,35, d. h. also dy = 10% 
= 1100 A. Es wurde nun Vo = 2,3. 
10% em! so gewählt, daß man sich 
in der Nähe des Punktes nd = 3/8) 
befindet, und die Konstante (* 
aus dem Anstieg der Kurve o () 
in der Umgebung von », ermittelt: 
C* = 1,15. Für die Schichtdicken- 
änderung Ad (in Ä) ergibt sich 


Abb.2. Die Schichtdickenverteilung einer 
Kryolith-Aufdampfschicht, gekennzeichnet 
durch Höhenlinien. Die eingetragenen Zahlen also nach (10), wenn man Aojo, 
geben die dort gemessene Dickenverminderung in % einsetzt 
in A gegeniiber der maximalen Schichtdicke 
(dy = 1100 A) an de 1100 _ de [A] 


Es wurde nun weiter auf vier, unter 45° stehenden Durchmessern im 
Abstand von 5 zu 5 mm das Reflexionsvermégen gemessen und Ad berechnet. 
In der Darstellung Abb. 2 wurde für die dickste Stelle der Schicht Ad will- 
kürlich gleich Null angesetzt und die Höhenlinien eingezeichnet. Aus dem 
einigermaßen glatten Verlauf der Höhenlinien — insbesondere auch in der 
Umgebung des Zentrums — und aus der Reproduzierbarkeit während der 
Messung kann man auf eine Genauigkeit von + 1 bis 2 A schließen, also eine 
Genauigkeit, die sonst nur durch Vielstrahlinterferenzen zu erreichen ist. 


4. Eine Anwendung zur Bestimmung der Verdampfungscharakteristik 


Der hier verwendete Verdampfungsofen bestand aus einem direkt geheizten, 
30 u dicken Wolfram-Blechstreifen. Er war so, wie in Abb. 3 angedeutet, aus- 
geschnitten und gebogen, daß eine einigermaßen gleichmäßige Aufheizung 
gewährleistet war. Der Restgasdruck während der Aufdampfung betrug etwa 
p = 3: 105 Torr (gemessen mit einem Ionisationsmanometer), so daß wir mit 
einer gradlinigen Dampfstrahlausbreitung rechnen können. 

Daß die dickste Stelle der Schicht exzentrisch liegt (Abb. 2), bedeutet, 
daß der Ofen keine zur Mittelnormalen der bedampften Fläche symmetrische 
Richtcharakteristik hatte. Weitere Versuche zeigten, daß diese Vorzugs- 
richtung der Öfen nicht gut reproduzierbar ist, sondern vermutlich von der 
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Lage der Substanz im Tiegel abhängt. . Trotzdem ergibt sich für den Abfall 

der Schichtdicke nach außen eine den theoretischen Formeln gut entsprechende 

Abhängigkeit. 
Für die Intensität des Molekülstromes, der auf 


einen senkrecht über dem Ofen im Abstand ry be- m ns 
findlichen, ebenen Auffänger unter dem Winkel 9 _ > 
auffällt (s. Abb. 3), gilt nach L. Holland) fiireine =| %| /r 
sogenannte punktförmige Dampfstrahlquelle 

und für eine flächenförmige Dampfquelle 


in 
(12b) 


J, ist hierbei die Intensität des Molekülstromes im 
Zentrum. Beachtet man, daß in unserem Fall r§ > a? 


ist 22cm; x = 0 = 3cm), so folgt für .die 
Punktquelle 


J “ae (=) 3a) — 


rel + (2/r9)*)"! Icm 


und für die Flichenquelle Abb. 3. Zur Charakteri- 
en stik des Verdampfungs- 
1 x \2 ofens 
(1 + (2/r9)*)° ro 

Die RENNER sind den Molekularstrahl-Intensitäten proportional; d.h. 

- =1-7, also (14) 
Der Anstieg 3/2 oder 2 der Kurve — “= = (= = liefert also die Entschei- 

0 


dung ob eine Punkt- oder Flächenquelle vorliegt. 

Aus der Abb. 2 wurde daher für verschiedene Richtungen der Durchmesser 
2x der in erster Näherung kreisförmigen Höhenlinien gemessen, gemittelt 
und — Ad/d, als Funktion von (z/r,)? 
in Abb. 4 aufgezeichnet. Es ergibt ad 


sich tatsächlich sehr gut angenähert 

eine Gerade, und zwar mit dem einen 

Anstieg 2,01 + 0,03. Es liegt also Fe 
ganz einwandfrei die Charakteristik Rn 
einer flächenförmigen Dampfquelle 
vor. Damit kommen wir zu dem 
gleichen Ergebnis wie L. Holland 0 0005 001 (£° 


und W. Stecklemacher'), die auf 
Eder Wege für M gF, „Schichten Abb. 4. Die relative Schichtdickenanderung 


als Funktion v ?.(22—= Durch 
und ähnliche Öfen ebenfalls Flächen- = 2: 


quellen fanden. Verdampfungsquelle/Schichtträger) 


4) L. Holland, Vacuum Deposition of Thin Films. London, 1956. 
5) L. Holland u. W. Stecklemacher, Vacuum London 2, 346 (1952). 
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Nat 5. Zur Anwendung des Verfahrens auf Mehrfachschichten 2 Schicht 
Für eine Zweifachschicht ergibt sich z. B. an Stelle von (3) u 
gehen 

R (®,; ®) = Ry + En. BAO, + 55.19; (15) [wie dur 

Ad, , 4v In d 

—=4n N] d, v AB, = —=4n N, (% Ad, + d;. Av) = = Einfachs 
d,. 0 

(16) 


Ay 
®,=Annd,v AD, = 47 Ng (v, Ad, + Av) = Do, ) 
Selbst wenn öR/ö®, und @R/2®, in einem gewissen Bereich um und | 
annähernd konstant sind, läßt sich das beschriebene Verfahren nur anwenden, 
wenn außerdem gilt P 


— (17) 


Es ergibt sich nämlich dann für » = », = const bzw. d = d, = const. 


18) also 


Hierbei 


an den Gin. (4): (5) entsprechende Dulane 6). Für 

Die Bedingung (17) bedeutet, daß die Schichtdickenänderung immer | 4 sind" 
proportional zur Gesamtdicke dieser Teilschicht sein muß. Derartige Schicht- 
dickenverteilungen werden aber praktisch häufig vorkommen, z. B. wenn die 
Verteilung durch die geometrische Anordnung von Verdampfungsquelle und | also in 
Schichtträger bedingt ist. In diesen Fällen ist das Verfahren natürlich nicht 
nur auf Zweifach-, sondern auch auf Mehrfachschichten anwendbar. 

Übrigens gibt es auch Mehrfachschichten, für die das Verfahren anwendbar Für 
bleibt, obwohl die Bedingung (17) nicht gilt. Das ist z.B. dann der Fall, 
wenn sich die Brechzahlen der Schichten nur wenig unterscheiden oder wenn 
die eine Teilschicht nur sehr dünn ist (n, d, <4/, A), wie man auf Grund von 
Näherungen zeigen kann. Für das Reflexionsvermögen schwach inhomogener 
Schichten schließlich gilt eine sehr ähnliche Formel wie für homogene Schich- 
ten’), in der nur die Fresnelkoeffizienten etwas anders definiert sind. Die 
Schichtdicke d und die Wellenzahl » gehen jedoch nur durch das Argument | !"4 da 
®—=4nndv in die Formel ein (R = mittlere Brechzahl). Es gelten also 
streng alle oben für eine homogene Einfachschicht abgeleiteten Beziehungen. 
— Kleine Störungen im Schichtaufbau, wie sie häufig vorliegen dürften®), 
haben also für unsere Betrachtungen keine Bedeutung. 


Da 


Nu 


6. Der Einfluß der Dispersion 


da aucl 
Bei der Ableitung der Formel (5) wurde vorausgesetzt, daß die Dispersion Bbung 
n (v) vernachlässigt werden kann. Das ist auch für reflexionsvermindernde PR 


6) Die Schichtdicken dı;o bzw. d2;9 können hier natürlich nicht mehr wie bei eirer 
Einfachschicht aus der Kurve R(v) bestimmt werden, sondern müssen anderweitig bekarnt 
sein, z. B. aus Messungen während der Herstellung der Mehrfachschicht. R 

7) H. Schröder, Ann. Physik (5) 39, 55 (1941). =e 

P. Bousquet, Ann. Physik Paris 2, 163 (1957). 


10) 
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Schichten recht gut erfüllt. Für reflexionserhöhende Schichten ergibt sich da- 
gegen zum Teil ein erheblicher Einfluß der Dispersion, der sich allerdings weit- 
gehend vermeiden läßt, wenn man Schichten der Dicke nd ~ 3/8 A benutzt, 
wie durch Betrachtungen — ähnlich?) — gezeigt werden kann: 

In der Umgebung der Wendepunkte ®, = 2/2; 5/22;... kann die Reflexion einer 
Einfachschicht durch folgende Formeln angenähert werden 

® = + AD cos (®, + = — sin ADR — AG 
ID = 4a(ndv — ny dy%) — 


ah 


R > 
ry „cos ® 


24 ,2 1 + 7? 


Hierbei sind A und B Funktionen von nity 
Für die Abhängigkeit R(Ad) (v und n konstant) folgt, da A und B unabhängig von 
d sind 


p also auch von » und damit auch von v. 


dR Ad 
= m =_— 2 
IR dd) |, Ad By Ad B, ®,y d,’ (21) 
also in Analogie zu (4) 

Für die Abhängigkeit R (Av) (d = const) folgt dagegen 

dR dA dB d(A®) a 


Da AD/, = O ist, verschwindet das mittlere Glied; weiter ist 


, dn „Av v, dn 
B, "dv is k Ip By (4 dy (m Av Ro C— (1 
und damit 


AR/Ry = (1 (25) 
CR, dv dv 
Nun ist entsprechend (20) are 
aA 2r, dern) 04 (ln) (26) 
er, + I Ory (1 + r 2 11? 


da auch bei reflexionserhöhenden Schichten r < 0,5 ist und damit 74< 1. Für die Um- 
gebung von ®, kann man wie in?) die Dispersion n(v) und auch die F Fresne Ikoeffizienten 


r(v) linear zu Ay ansetzen und mit (26) dA/dv ausrechnen: By vr Elan 
wor dr +2r dry 2 


®) G. Koppelmann u. K. Krebs, Z. Physik 145, 486 (1956). 
) Das Zeichen /,deutet an, daß der betreffende Wert an der Stelle v,; dyzu nehmen ist. 
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Für Dy = an; 7/.2;... gelten die gleichen Formeln; nur besitzt B und damit auch 
C das umgekehrte Vorzeichen, da cos (®, + J®) + sin A® in (20) einzusetzen ist. 
Zusammenfassend läßt sich also aussagen: Im Falle nicht vernachlässigbarer Di- 

spersion gelten statt (4) und (5) die Formeln (22) und (25) 

AR Ad AR ‚Av . -C 1 dA dn 

— = — = C’—,_ wobei ; — % +4 

Ry dy Ry % ( C Ry dv Ny dv" 
Es ergibt sich also fiir relative Differenz (C’ — C)/C die Summe von zwei Termen. Der 
erste rührt daher, daß sich mit n (v) die Reflexionskoeffizienten r,; 7,,; andern. Der zweite 


(28) 


dn 
Term beruht auf der Tatsache, daß in AD = 42d, (m + — ») A» die Disperison eben- 


dv 
falls eingeht. Je nach der Dicke der Schicht n d = 1/, A; A; P/g A; .. . baw. A: 7/g Ar... 
verstärken oder schwächen sich die beiden Anteile. 
Als Beispiele wurden die Größen (C” — C)/C mit den Konstanten einer ZnS-Schicht 
(n = 2,39 dn/dv = 0,21 - 10-4 em) und einer MgF,-Schicht!!) = - 1,40; dn/dv = 0,015-10-4) 


auf Kronglas (n = 1,52; dn/dv = 0,018 - 10%) berechnet (% = 2-104 cm): 

ZnS (-0)/C 0,59 0,04 0,26 0,11 0,22 0,14 
MgF, —0,14 +0,07 +0,003 +0,03 


Während also die relativen Differenzen im Falle der MgF,-Schichten (außer für 
nd = !/, 4) praktisch bedeutungslos sind, erweisen sie sich für die ZnS-Schichten als z. T. 
sehr hoch. Für ?/, A kompensieren sich dagegen die Einflüsse fast vollständig. Man kann 
zeigen, daß diese Tatsache nicht auf das spezielle Zahlenbeispiel beschränkt ist, sondern 
für stark reflexionserhöhende Schichten unabhängig von derGröße der Dispersion allgemein 
gilt. Es ist also möglich, auch bei reflexionserhöhenden Schichten aus ziemlich stark di- 
spergierenden Substanzen das Verfahren der Messung der Schichtdickenverteilung ohne 
Korrekturen anzuwenden, indem man %/, A-Schichten benutzt. 


Herrn Prof. K. Krebs danke ich für zahlreiche Diskussionen und wert- 
volle Ratschläge. Der Deutschen Forschungsgemeinschaft bin ich für finan- 
zielle Unterstützung dieser und weiterer Arbeiten zu besonderem Dank ver- 
pflichtet. 


1) Als Beispiel für eine reflexionsvermindernde Schicht wurde eine MgF,-Schicht 
gewählt, da für Kryolith in der Literatur keine Angaben über die Größe der Dispersion zu 
finden sind. Die Verhältnisse werden jedoch hier — wie auch bei den anderen niedrig- 
brechenden Substanzen — sehr ähnlich liegen wie beim MgF,. 


he Berlin, II. Physikalisches Institut der Technischen Universität Berlin. 


maya Bei der Redaktion eingegangen am 31. Juli 1959. Bellen 4 
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Fourier-Analyse des Mikrofeldes 
in einem Plasma. II 


Inhaltsübersicht 


In einer ersten Mitteilung!) wurde eine Fourier- Analyse des zeitlichen 
Verlaufs einer beliebigen Komponente des elektrischen Mikrofeldes in einem 
Plasma durchgeführt. Dabei wurde von der Coulomb-Wechselwirkung der 
geladenen Teilchen des Plasmas abgesehen, ihre Bewegung also als gleichförmig 
und geradlinig betrachtet. Da diese Annahme bei größerer Dichte der geladenen 
Teilchen nicht mehr zulässig ist, wird hier die Wirkung der stark gekrümmten 
Teile der Elektronenbahnen auf das Frequenzspektrum untersucht. Es ergibt 
sich, daß diese Bahnteile hauptsächlich Beiträge zu den höchsten Frequenzen 
liefern. Für die Wahrscheinlichkeit der Fourier- Amplituden der Feldstärke 
ergibt sich auch hier eine Gauß- Verteilung, die für hohe Frequenzen gegen- 
über der in I gefundenen nach größeren Amplituden verschoben ist. 


Die Fourier- Analyse einer Komponente des elektrischen Feldes in einem 
Plasma wurde in I für den Fall durchgeführt, daß die Bewegung der geladenen 
Teilchen als geradlinig und gleichförmig angesehen werden kann. Hier soll 
der Einfluß der Ionen auf die Elektronenbewegung berücksichtigt werden. 

Wir betrachten zunächst die Hyperbel- 
bewegung eines Elektrons (Ladung — q) unter 
der Wirkung der Coulomb-Kraft eines als 4 
ruhend gedachten Ions (Ladung Z q). Legt man 
den Anfangspunkt eines Koordinatensystems Z, 7 
in den Asymptotenschnittpunkt und nimmt die 
Hyperbelachse als %-Achse, so läßt sich die Be- 
wegung mit Hilfe eines Parameters £ durch 


darstellen. Hierin ist a = as die Halbachse, ¢ die 


Exzentrizität, r = a ®—1 der StoBparameter, 


') G. Hettner und H. Wagner, Ann. Physik (7) 
4, 89 (1959), im folgenden als I zitiert. In der letzten Abb. 1. Hyperbelbahn eines 
Formel dieser Mitteilung ist leider ein Fehler stehen- Ejektrons. I Ion, a große 
Halbachse, & Exzentrizität, 


geblieben, es muß heißen: 0 für w 
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v die aa Geschwindigkeit des Elektrons (vgl. Abb. 1) und 
v_ mv 


m =2=F5. () 


Um # und g Fourier-integrabel zu machen, fügen wir den rang ” ei 
(x reell > 0) hinzu und erhalten = 
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und mittels partieller Integration 
a{ 2a (-) i 
+ 
Tu +0? ' w-ia 
In ähnlicher Weise ergibt sich 


IO Tor 
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co 
- fem dé und I” = fe 
6 


ein, so findet man 


(-) Wy 1 (+) i Wo 
= — d 
Is e(w — 124) € und e(m +10) 
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4 Vel. G. Hettner, 2. 150, 
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7 > % 4 


Für x > 0 wird > 


+00 .@ 7 

mr 
+ls -le=— < Sin ge ® 
also 
— , amd ar — ri 


Da geradlinige Elektronenbahnen in I behandelt wurden, betrachten wir hier 
nur die stark gekriimmten Teile der Hyperbelbahnen, für die | | << 1 ist. Wenn 
die Beträge der Integranden in J~ und J“) für die nahezu geradlinigen Teile 


der Bahn klein sein sollen, muß en 


sein. Diese Bedingung ist erfüllt für 


Die Grenzwerte (9) werden also nur dann für uns brauchbar sein, wenn die 
Annahmen « = 0 und x = a, zu nahezu denselben Werten für die Amplituden _ 
und führen. Das ist für 


der Fall. Diejenigen Teile der Integrale IC) I) (3), die: sich Inter- 
vall 1 << |&| < co beziehen, gehen nämlich sowohl für « = a, wie fiir a = 0 


0 
Xp 


0) 
mit wachsendem — gegen z 
0 


einem Lemma von Riemann- 


Lebesgue?°)), und auch die übrigen th Wr 


in (6) und (7) enthaltenen Ausdrücke 
nehmen in beiden Fällen die gleichen P 
Grenzwerte an. Für werden 

dagegen die Amplituden und yo at 
proportional @ und verschwinden 
daher gegen die Beiträge der geraden 

Bahnen, die für kleine @ vonwun- 


abhängig werden. Wir verzichten 
daher auf eine Berechnung der Bei- > u 
träge der Bahnkrümmungen zu den ert) 


Amplituden kleiner Frequenzen. » aise 


Um überdie Elektronenbahnen u 
summieren, führen wir 2 kartesische Abb. 2. Koordinatensysteme 2 yz und 
x’ y’z’. I Ion, P Asymptotenschnittpunkt, 


®») Vgl. z. B. Whittaker-Watson, x Azimuth der Bahnebene gegen die Meri- — 
Modern Analysis, Cambridge 1952, S. 172. dianebene P z. 


| 


(4) 
0 = € 
dé 
(5) 
"4 
(7) 
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Koordinatensysteme x,y,z und x’, y’,z’ mit einem Ion J als Anfangspunkt und bi 
ein (Abb. 2), deren gegenseitige Lage durch die Winkel # und gegeben sei. 
Die z’-Achse sei die Hyperbelachse einer Elektronenbahn, ihre Ebene habe q 


gegen die x’ z’-Ebene das Azimuth x. P sei der Schnittpunkt der Asymptoten, 
also der Anfangspunkt der Koordinaten 7, 7 (Abb. 1). Die Komponenten der 
Fourier-Amplituden in den Achsenrichtungen dieser Koordinatensysteme | Aus de 
werden wegen 2’ = cosa, y = 7 sina, =ea—F und w> 0 


Io = Yo COSA, Yo = YoSNA, = — Low (12) 
und = Bezeic! 
Be = cos — Yo SiN + COS sin d 
| — 0 (13 


Da der Perihelabstand des Elektrons vom Ion im allgemeinen klein gegen den 
Abstand des Ions vom Aufpunkt ist, beschreiben wir die Wirkung des Elektrons 
im Aufpunkt als die Wirkung von 2 Oszillatoren mit den Amplituden x, und 
2», die im Ion lokalisiert sind. Ist r, der Abstand des /-ten Ions vom Aufpunkt, 
so werden die Komponenten der Amplitude der elektrischen Feldstärke im fda das 


2 
= ley = 0, Jone hes (14) 
fr 


Wir führen nun eine beliebige, aber feste Z-Achse mit dem Aufpunkt A als 
Anfangspunkt ein (Abb. 3). Die z-Achse bilde mit ihr den Winkel #,, die I 
x-Achse liege in der zZ-Ebene. Dann wird die Z-Komponente der Feldstärken- 
amplitude, hervorgerufen durch das k-te Elektron solange 
wi 2 unter dem Einfluß des /-ten Ions 
fu" = fo,z 6089, + foo,x sin 

und wir setzen 


ak: — Re = In 16) 


_— k-te Elektron passiere zur Zeit ¢ = t, das Perihel, | 
fa ist (vgl. I) 


A v(t 19 ) ei ot dt 
Abb.3. A Aufpunkt, 22 _ und da 
IIon, f,,f, Komponen- iwt, +0 
ten der Feldstärke in f a(t) = 0) 
Elektron unter der : 
Wiens und ebenso > vat Da die 
Ga 91» Og al 

Zur Berechnung der Wahrscheinlichkeitsverteilung der Feldstärkenamplituden die je 
verwenden wir wie in I die Methode von A. Markoff. Es seien p, und p, On 


die Gesamtheit der geeignet transformierten Lagen- und Impulskoordinaten 
des k-ten Elektrons und des l-ten Ions (k = 1, 2,..., N’; 1=1,2,...,N”), | 
t(p,) dp, und o(p,) dp, die Wahrscheinlichkeiten für die Intervalle p, ---p,-+ [Far die 
dp, und p,---~, + dp,. Wir setzen tes Ons 9 

l naten r 


| 


le 


Up I 
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kt [und bilden das Integral 

a] pin 

er 4 (N’’) 1 


ne | Aus denselben Gründen wie in I können wir die e-Funktion durch cos u*:! 


und den Cosinus durch 1— —-(u*!)® ersetzen. Daher wird, wenn wir die 


2 
Bezeichnungen 
einführen, 
nd o(p,) dp, 


m | da das Integral über das Produktglied wie in I verschwindet. 
Die Mittelung über die Koordinaten p,, d. h. über alle Ionenlagen, ergibt 


4) (22 
I o(pı) dp, 1 v” (ai v2) 
us 44 —— 92 + — 1 9? + 
on [solange N” (---)< 1. Mit den Bezeichnungen 
5) 2 = ' (23) 
(k) (k) 
wird 
6) Nr 
(é) N 2 2 
— ay oT + bu 05 (24) 
| t(p.) dp, ( 
(k) 
und damit schlieBlich 
7) 
Ty’ x” (01: 02) = (a2, + be 03), ). (25) 
Da die e-Funktion für o,, 0, > co gegen 0 geht, dürfen wir (25) für beliebige 
9,0, als gültig annehmen. 
a Wir betrachten die Ionen als ruhend und berücksichtigen nur diejenigen, 
p die sich in der Debye-Kugel um den Aufpunkt befinden. Dann wird 
o(p,) dp, = r? dr, sin 9, dd, dq, J o(p,) dp, (26) 


+ #Für die Elektronen benutzen wir dieselben Koordinaten wie in I, nämlich 
Tes Dn, Pps &, |, mit dem Ion als Anfangspunkt. Allerdings sind diese Koordi- 
(8) Jnaten nur für die nahezu geradlinigen und gleichférmig durchlaufenen Teile 
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der Hyperbeläste verwendbar, weil nur für diese der Elektronenort allein 
durch die Koordinate /, bei konstanten übrigen Koordinaten beschrieben wird. 
Da r, in diesem Falle das Lot vom Ion auf die Asymptote ist, ist r, identisch 
mit dem hier unter derselben Bezeichnung eingeführten Stoßparameter. Wir 
können diese Koordinaten trotz ihrer Einschränkung auf bestimmte Gebiete 
des Ortsraumes auch hier verwenden. Da nämlich D stets wesentlich größer 
ist als a, können wir eine Länge B definieren, so daßa < B< D. Dann liegen 
innerhalb der Kugel vom Radius B mit dem Ion als Mittelpunkt die stark 
gekrümmten Teile der Elektronenbahnen. Wir schließen nun bei der Inte- 
gration über die Ortskoordinaten der Elektronen zur Zeit t = 0 diese Kugel 
aus. Das bedeutet, daß Elektronen, die in einem Zeitintervall von der Größen- 


ordnung — 2< t< + 2 den stark gekriimmten Teil ihrer Bahn durch- 
laufen, nicht mitgezählt werden. Da im ganzen alle Elektronen berücksichtigt 
werden, die im Zeitintervall —T <t<+T mit 7 > [vgl. 1(10)] die 


Debye-Kugel passieren, führt die Nichtberücksichtigung der Kugel vom 
Radius B zu keinem wesentlichen Fehler. Wir dürfen dann auch wie in I 


i 


und [I (15)] 
cos "di, = f dl, = R, [ cos* „sin dl, = 0 
- Lk -Ix 
setzen. 


Da wir die asymptotischen Geschwindigkeiten v und damit die Gesamt- 
energien aller Elektronen einander gleich annehmen, so ist der Boltzmann- 
Faktor eine Konstante und hebt sich bei der Mittelbildung ebenso wie die 
Masse und die Geschwindigkeit heraus. Wir können daher ähnlich wie in I (2) 
setzen 

t(p,) dp, = r, dr, sin 9, dd, dl,dp,da,; | t(p,) dp, = 82? B®R. (29) 
Um nun die in (25) auftretenden Größen a?, und 2, zu berechnen, drücken 
wir zunächst die Quadrate der in (16) definierten Größen at! und ph mit 
Hilfe von (15), (14), (13), (12), (17) und (9) durch die Integrale /g und Ig 
und durch die Ionen- und Elektronenkoordinaten aus: 


2 

= [ cos? {cos = alg-2cosd, + sin r,, Ig cos 2 sind 

cos sin {- - -} + sin? d, Icos "als (— cos p, sin 9,) (30) 

2 

+ sin = r,, Ig (cos x, cos p, cos #,, — sin x, sin 2 ] ‘ 
Bezeichnen wir mit s’ den kleinsten Wert von r,, bei dem die Darstellung des 
Feldes der Elektronen durch im Ion lokalisierte Dipole [Gl. (14)] noch sinn- 
voll ist, so ergibt die Integration über r, 
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Wegen B < D ist es hierbei unwesentlich, daß die Kugeln mit dem Radius B 
für die äußersten Ionen nicht voll in das Integrationsgebiet fallen. Daher 
wird nach (20), (26) und (30) und nach Integration über 3, und 9, wen 


(ky? _ 2 1% cos?” £5 eos? 2 2 
cos? #,, + 
= 08 os? #,, + sin? 3, cos? p,} (32) 


+ rf 1% sin?® sin?d, cos? &, +- cos? 9, cos? &, cos* + sin? x, sin?Q, 
ol, . ol, 
— 2 cos 9, cos &, sin &, Cos sin y,} + ar, Ig cos 3] 


In pi” ist sin = mit cos = vertauscht. Nach (23) und (29) folgt hieraus 
schließlich durch Integration über J, [Gl. (28)] und über 9, und ¢, 


= ba (a? I& + rt Ic) r.dr,. (33) 
Hierin sind 7, und r, der kleinste und ne größte in Betracht kommende Stoß- 
parameter. Wir bestimmen ihre Werte dadurch, daß wir den kleinsten Perihel- 
abstand gleich dem Radius s der Ionen, den größten gleich B setzen. 

In (33) führen wir als Integrationsvariable die Exzentrizität & ein, setzen 
also rf = (e® — 1) a? und r, dr, = a? ede und erhalten 

Ya 
(a? 1% + rf 1%) r,, dr, = at f (1% + (e — 1) 18 ede. (34) 


& und & sind durch (e,— 1)a=s und (e,— 1)a= B bestimmt, also ist 
14+ g=1+2>1. (35) 


Um die Integration in (34) auszuführen, benutzen wir die für eine beliebige 
Zylinderfunktion Z,(z) gültige Beziehung 


2, Zp) + (36) 


die unmittelbar aus der Besselschen Differentialgleichung folgt. Setzt man 


in (36) Z,(2) = HY (pe), so wird aot 
€) + (. 1) Hy” (pe)\ ede = = peH, (pe) HS (pe) 
oder für p = auf Grund der Beziehungen H (17) ini iy. wi 
1) Ig} ede = Isle (37) 
oder nach (34) und (35) ~~, 
(a? 1%; + IR) dr, = Is «(2 | Gk 
(38) 
= gf (w, &) — y(o, &)), 
worin 
y(o,e) = ls (=..). (39) 
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Wegen & > 1 ist y(o,&) in (38) zu vernachlässigen. (33) und (38) ergeben 
schließlich 

gaa 

a, = ba = 622 D BY &)- (40) 
Sind », und n, die Anzahlen der Elektronen und Ionen in der Volumeneinheit, 
so ist die Zahl der Elektronen, die im Intervall — T <t<-- T die Kugel 
vom Radius B um ein Ion passieren, nach I (12) 


(41) 
und die Zahl der (ruhenden) Ionen in der Debye-Kugel um den Aufpunkt 


4x 


also 


2 2 be = &ı) (43) 
unabhängig von D und B; folglich ist nach (25) zn i bain, DR 


Ty (0p 09) = 
Ebenso wie in I (22) und (23) ergibt sich die Wahrscheinlichkeit, «laß der 
Betrag |F „| der Feldstärken-Amplitude zwischen | Fo und | po 
liegt, 


re 


| wie ie Der wahrscheinlichste Wert der Feld- 

&=1 stärken-Amplitude ist also 2%. Eine 


N ug Reihe von Werten der Funktion y(o, &) 
S wurden nach (39) und auf Grund der 
Darstellung der Integrale /g und Jz in H 


(16) mit der PERM berechnet) (© 
0,5---10; e=1--- 10). In Abb. + ist 


ein Teil der Ergebnisse aufgetragen, 


0 vone nehmen die Funktionswerte weiter 
o stark ab. Man erkennt, daß die Elek- 
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(39) definiert. Für größere Exzentri- 
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weitem die größten Beiträge liefern. Für diese ist der Perihelabstand 
(e—1)a<a. Man kann daher in (46) s’ =a setzen und erhält 
2n v \3 p(a, &) 
4 ‘ 


e Für große Werte von — und ¢ = | läßt sich ein Tales ert für y ableiten: 
% 
4nu t 
für 1. (48) 
1) V3 @ 


kt | Die mit der PERM berechneten Werte zeigten, daß dieser Näherungswert 
[477 > - . . 
für —— 1 um etwa 16% und für — = 5 um etwa 6,5% zu klein ist. 

Wo 0 o /o 


0 
Da g(w) und damit auch @ für große w exponentiell gegen 0 geht, y(w, 1) 


aber wie w~! (48), Y also wie w~* (43), so gibt es ein Gebiet hoher Frequenzen, 

3) B in dem Y>@ ist. Auf diesen Teil des 
Spektrums haben also die Bahnkriim- wey 
mungen einen wesentlichen Einfluß. Abb. 5 

4) zeigt qualitativ die Wahrscheinlichkeits- 

verteilungen Wy und Wr für ein gentigend I 

0, | großes wm. Wy bezieht sich auf die in I gt 

behandelte gleichförmig geradlinige Be- Wy 
wegung der Elektronen, Wy, auf die hier 
untersuchten gekriimmten Teile der Bahnen. 

15) Die durch die beiden Teile der Bahnen Abb. 5. Wahrscheinlichkeitsvertei- 
hervorgerufenen Feldstärken-Amplituden lungen für die Feldstärken-Ampli- 
FR=Ar+iBı und Fır= Ari Bir  tude|F füreine hohe Frequenz o. 

sind als statistisch unabhängige Ereignisse MW, gilt für die nahezu g geraden Teile, 

!6) | anzusehen, weil die Bahnkrümmungen von W,, für die stark gekrümmten Teile 

Id. den statistisch verteilten Lagen der Ionen der Elektronenbahnen (qualitativ) 


abhängen, die geraden Teile aber nicht. 

Man erhält daher die Wahrscheinlichkeitsverteilung W (A, B) für A = Aı + Ary 
€) | und B= Bı + By durch eine Faltung aus W(A,, By) und W(Aır, By). Aus 
den beiden Gauß- Verteilungen entsteht wieder eine Gauß-Verteilung. Für 


H | den Betrag |F,„ der resultierenden Feldstärkenamplitude F, erhält man die 
Verteilung 
(0) |2 


2(®+%) 


© München, Institut für Theoretische Physik der Technischen Hochschule. = 
ter Bei der Redaktion eingegangen am 11. August 1959. 
k gegang g 
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der 
aus der positiven Säule von Sauerstoff-Glimmentladungen 
effundierenden und extrahierten Ladungsträger 


Von Günter Brederlow 

Inhaltsiibersicht 


Mit einem einfach-fokussierenden Massenspektrometer geringen Auf- 
lösungsvermögens aber hoher Empfindlichkeit werden die aus der stationären 
positiven Säule von Sauerstoff-Glimmentladungen effundierenden und extra- 
hierten positiv und negativ geladenen Teilchen in der T- und H-Form untersucht. 

Die bei Entladungsstromstärken von 20 bis 120 mA und einem Entladungs- 
druck von 2,5 Torr gemessenen Teilchenströme werden mitgeteilt und ver- 
glichen. Das Verhältnis der Ströme ist von der Entladungsform abhängig. 

Die Größe des Effusionsstromes wird wesentlich durch die Amplitude der 
Potentialschwankungen des Plasmas vor der Effusionsöffnung bestimmt. Die 
Potentialschwankungen treten sowohl als statistische Schwankungen als auch 
als Schwingungen auf. 

Mit einer Gegenfeldmethode wird die Energie der effundierenden Teilchen 
gemessen. Es kann bei Potentialschwingungen ohne überlagerte statistische 
Schwankungen gezeigt werden, daß Proportionalität zwischen der Amplitude 
der Potentialschwingungen und der Teilchenenergie besteht. Aus den MeB- 
ergebnissen folgt, daß die Schwingungen vor der Effusionsöffnung mit einer 


Komponente senkrecht zur Rohrwand auftreten. Re u 
Im. 


= 


In Sauerstoff treten zwei Entladungsformen auf, die sich in wesentlichen 
Säulenmerkmalen unterscheiden. Diese Erscheinung wurde von Seeliger und 
Mitarbeitern!) untersucht und wegen der charakteristischen Differenz in der 
Höhe des Gradienten mit H- und T-Form bezeichnet. Die jeweilige Ent- 
ladungsform brennt nur in einem bestimmten von der Rohrweite abhängigen 
Stromstärke- und Druckbereich stabil. Die H-Form bei geringem Druck und 
hohen Stromdichten, die T-Form gerade umgekehrt. Der Übergang erfolgt 
sprunghaft. Das Übergangsgebiet verschiebt sich mit abnehmendem Druck 
zu kleineren Stromdichtewerten, wobei bei zu tiefem Druck schließlich nur 
die H-Form auftritt tnd bei zu hohem entsprechend nur die T-Form. 

In der vorliegenden Arbeit soll versucht werden, nach einem von Pahl?) 
angegebenen Meßverfahren durch massenspektrometrische Untersuchung der 


1) R. Seeliger u. A. Wichmann, Ann. Physik 9, 235 (1951); R. Seeliger, Z. Natur- 

j forschg. 8a, 74 (1953). 
2 7 N 2: 6: y 57 

“te ) M. Pahl, Z. Naturforschg. 12a, 632 (1957) 
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durch eine kleine Öffnung in der Rohrwand austretenden positiven und nega- 
tiven Träger einen Beitrag zur Klärung dieser Erscheinung zu geben. 

Die effundierten bzw. abgesaugten Teilchen werden beschleunigt und in 
einem für diese Zwecke entwickelten stigmatisch fokussierenden Massen- 
spektrometer geringen Auflösungsvermögens aber hoher Empfindlichkeit 
nachgewiesen. 

Das Entladungsrohr 


Die Entladung brennt in einem Glasrohr, dessen Aufbau und Abmessungen 
aus der Skizze (Abb. 1) ersichtlich sind. Als Kathode dient ein Nickelbecher und 
als Anode eine Nickelscheibe. Der mit der Effusionsöffnung versehene Teil 
des Entladungsrohres befindet sich in einer evakuierbaren Kammer. Die 
Verbindung erfolgt durch einen mit Silicon-Hochvakuumfett gefetteten 
Schliff. Diese Abdichtung 
erweist sich auch bei höch- 
ster thermischer Belastung 
als sehr zuverlässig und \ 
sicher. 

Die Effusionsöffnung im 
Entladungsrohr wurde 
durch Atzen mit Flußsäure kee 
mittels einer feinen Nadel 
in einem möglichst dünn- 
wandigen Teil der Wand Bihar 
hergestellt; sie war stets 
triehterförmig. Andereng- 
sten Stelle im Innern der 9 2 46m 
Rohrwand hatten die ver- 
schiedenen Öffnungen einen 
Durchmesser im Bereich *bb.1. Skizze des Entladungsrohres und des Massen- 

4 spektrometers ohne Absaugfolie. A Ablenkkammer, 
von 60—100 u. Sollten ge- 


zn B Beschleunigungskammer, C Auffänger, D Effusions- 
ladene Teilchen abgesaugt öffnung, M Magnet 


werden, dann wurde statt 

dessen am Rohr eine Platinfolie von 20 Dicke aufgeschmolzen. Diese 
ist mit einer Öffnung von 58 u Durchmesser versehen und bildet eine ungefähr 
0,3 mm? große Fläche der inneren Rohrwand. Zwei weitere im gleichen Ab- 
stand von der Anode aufgeschmolzene Platinfolien mit einer Gesamtfläche von 
imm? liefern das Bezugspotential für die Absaugspannung und eine vierte 
das für die Beschleunigung der austretenden geladenen Teilchen. 

Durch eine Vorrichtung zur Stabilisierung des Druckes, über die an anderer 
Stelle ausführlich berichtet wird®), wird der Druck im Entladungsrohr nahezu 
konstant gehalten. Die Druckschwankungen im Entladungsrohr sind kleiner 
ils 0,1 Torr. Wegen der Schliffverbindungen konnte das Entladungsrohr nicht 
ausgeheizt werden. Das Rohr wurde aber mit Strömen bis zu 300 mA unter 
häufigem Gaswechsel eingebrannt. Der Sauerstoff wurde durch Erhitzen 
von reinstem Kaliumpermanganat gewonnen. 

An die Beschleunigungskammer schließt sich das Massenspektrometer an. 
Die Verbindung erfolgt über einen Planschliff mit zwischenliegendem Dich- 


G. Brederlow, Vakuum-Technik 4, 108 (1959). 
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tungsring aus Gummi. Das in die Beschleunigungskammer durch die Effu- 
sionsöffnung strémende Gas wird mit einer Hg-Diffusionspumpe abgesaugk 
Der Enddruck ist stets kleiner als 10 Torr. 


Der Aufbau des Massenspektrometers ist der Abb. 1 zu entnehmen. Da 
bei der Effusion der negativen Ionen nur sehr geringe Teilchenstréme zu er- 
warten waren, mußte ein Gerät mit hoher Empfindlichkeit unter Verzicht auf 
großes Auflösungsvermögen entwickelt werden. Beide Eigenschaften sind bei 
dem vorliegenden Verwendungszweck vereinbar; denn schon mit einem Auf- 
lösungsvermögen von 4 können eventuell vorhandene Ionen der Masse 64 
und 48 getrennt werden. 

Eine Arbeit von Schönheit?) gab die Anregung zum Bau eines einfach- 
richtungsfokussierenden Massenspektrometers, das im wesentlichen aus einer 
elektrostatischen Rohrlinse, dem Magneten und dem Auffängersystem besteht. 
Die nachzuweisenden geladenen Teilchen fokussiert und beschleunigt eine 
elektrostatische Rohrlinse. Das Radienverhältnis und die Abmessungen der 
Rohre sind so gewählt, daß das der Effusionsöffnung näherliegende Rohr ein 
Potential erhält, das in seiner Polarität dem des anderen entgegengesetzt 
gerichtet ist. 

Das Massenspektrum wird durch Veränderung des Magnetfeldes abge- 
tastet. Um ein möglichst homogenes Magnetfeld mit geringem Streufeld zu 
erhalten, sind zwei Eisenplatten in die Wandungen des Massenspektrometers 
eingelassen worden, so daß ein Polabstand von 6 mm erreicht wird. Die Ab- 
dichtung erfolgt durch Apiezonwachs. Um Störungen der Ionenbahnen in 
der Beschleunigungskammer durch das Streufeld des Magneten zu vermeiden, 
wurden die Teile der Rohrlinse aus Eisen angefertigt und eine zwischen- 
liegende Eisenplatte angebracht. Bei dieser Anordnung konnte durch eine 
empfindliche Hallsonde kein Streufeld nachgewiesen werden. Die elektrische 
Durchführung zur Rohrlinse erfolgt durch eine Schliffverbindung. 

Der Ionenstrom kann sowohl in seiner Gesamtheit in einem in der Spektro- 
meterachse befindlichen Auffänger als auch in seinen Ionenanteilen in einem 
20 mm von der Achse entfernt liegenden Auffänger gemessen werden. In den 
ersteren gelangen die Ionen durch eine Öffnung von 10 mm Durchmesser und 
in den zweiten durch einen Spalt von 1,5 mm Breite. Durch diese Anordnung 
kann kontrolliert werden, ob im Massenspektrum alle Ionenanteile im vollen 
Betrag erfaßt werden. Beide Auffänger sind mit einer Elektrode, die durch 
ein negatives Potential das Entweichen der Sekundärelektroden aus den Auf- 
fängern verhindern soll, in eine Glasglocke eingeschmolzen, die über einen 
Planschliff mit einem Dichtungsring aus Gummi am Massenspektrometer be- 
festigt wird. Um den Auffänger gut zu isolieren, wurden möglichst lange 
Glaswege gewählt. Der Auffänger kann gegen einen Leuchtschirm ausgewechselt 
werden, der es ermöglicht, den Ionenstrahl visuell zu beobachten. 

Durch die Potentialschwankungen des Plasmas bedingt, haben die aus- 
tretenden geladenen Teilchen eine erhebliche Energiebreite. Es war daher 
nicht möglich, den Ionenstrahl so zu fokussieren, daß sein Durchmesser im 
Austrittsspalt kleiner als 1,3mm war. Bei einer Beschleunigungsspannung 


4) E. Schönheit, Z. angew. Physik 9, 171 (1957). 
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von 1800 Volt konnte daher nur ein Auflösungsvermögen von 6 erreicht 
werden. Das Massenspektrometer wurde aus einem Aluminiumblock herge- 
stellt. Seine Ablenkkammer wird mit einer Hg-Diffusionspumpe evakuiert. 
Der Druck wird mit einem Ionisationsmanometer gemessen. Sein Wert ist 
stets kleiner als 10° Torr. 

Der elektrische Aufbau der MeBanordnung = = ~~ 


Beim Entwurf der elektrischen Schaltung der Meßanordnung mußte beach- 


tet werden, daß das Massenspektrometer zu erden ist und für die Beschleuni- 
gung der geladenen Teilchen ein möglichst festes und zweckmäßiges Bezugs- 
potential geschaffen wird. Die aus Abb. 2 ersichtliche Schaltung erfüllt diese 
Forderungen. Die einfachste Lösung wäre, an der der Effusionsöffnung 
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Abb. 2. Schaltskizze der Meßanordnung mit Absaugfolie. A Absaugfolie, B Folien fiir das 
Bezugspotential, € Rohrlinse, D Auffänger, KB Kompensationsmeßbrücke 


gegenüberliegenden Folie, die das Wandpotential der Umgebung annimmt, 
die Beschleunigungsspannung zu legen und so die gesamte Meßanordnung auf 
das entsprechende Potential zu heben. Leider ist dies nicht möglich, da über 
die Isolationswiderstände in dem Entladungsstromkreis, die in der Skizze 
summiert durch einen gestrichelten Widerstand eingezeichnet sind, ein Strom 
zur Erde abfließt. Dieser Strom, der auch über die Platinfolie in die Entladung 
oder umgekehrt fließen muß, bewirkt aber, daß die Folie einen Potentialwert 
annimmt, der von dem der Wand verschieden und daher als Bezugspotential 
ungeeignet ist. Durch die aus der Schaltskizze ersichtliche Änderung kann 
dieser Mangel beseitigt werden. Die dem Entladungsrohr parallel geschaltete 
Triode und der Widerstand von 2 MQ bilden einen Kathodenverstärker, dessen 
Gitter über einen Schutzwiderstand mit der Folie verbunden ist. Da Gitter 
und Kathode eines solchen Verstärkers stets nahezu den gleichen Potential- 
wert annehmen, liefert das Kathodenpotential der Triode einen Bezugspunkt 
für das Beschleunigungspotential, ohne daß ein nennenswerter Strom über 
die Folie abfließt. Die Potentialschwankung des Plasmas, die in Abhängigkeit 
von der Entladungsstromstärke sowohl als statistische Schwankungen als 
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auch Schwingungen mit einer Frequenz bis zu 40 kHz auftreten, bewirken, 


daß das Potential der Folie und dementsprechend auch das der Kathode der we 
Triode in ungiinstigsten Fallen um einen Betrag von 40 Volt schwankt. Das oo 
feste Potential des Massenspektrometers, das auch das Potential des einen D 
Rohres der Rohrlinse ist, erfordert aber zwecks gleichbleibender Fokussierung El 
eine feste Spannungsdifferenz zwischen beiden Rohren und demzufolge auch | Weiss 


ein festes Bezugspotential. Durch ein RC-Glied, das parallel zum Kathoden- 
widerstand liegt, wird ein festes Potential gegeben. Es entspricht praktisch 
dem Mittelwert des Wandpotentials in der Umgebung der Effusionsöffnung 
bzw. der Absaugfolie und eignet sich daher als Bezugspotential. Es besteht D 
jetzt zwischen Rohrlinse und Effusionsöffnung oder Absaugfolie keine feste | Schw 
Potentialdifferenz. Die nachzuweisenden Teilchen haben, wie auch Gegenfeld- | breite 
messung von Weimer®) in Edelgasentladungen ergeben haben, eine Energie- | missi 
breite, die nahezu den Schwankungen des Plasmapotentials entspricht. Im D 
vorliegenden Fall dürfte dieser Wert maximal 40 eV erreichen. Diese Energie. | |, ng 
inhomogenität beeinträchtigt die Fokussierung des Ionenstrahles und damit gezeig 
auch das Auflösungsvermögen des Massenspektrometers. ER 

Die Potentialschwankungen der Folie bewirken, daß ein Wechselstrom | deute 
von der zwischen Gitter und Folie liegenden Verbindungsleitung über die | mach 
Leitungskapazität abfließt. Durch eine von Attree®) angegebene Ab-f 0,1 T: 


schirmung, die in der Skizze gestrichelt eingezeichnet ist, konnte diese störende | haber 
Kapazität fast vollkommen beseitigt werden. mitte 

Sollen geladene Teilchen abgesaugt werden, so wird statt der Effusions- | wird. 
öffnung in der Glaswand eine durchbohrte Platinfolie mit möglichst kleiner Ai 


Wandfläche aufgeschmolzen. Das Bezugspotential für die Potentialänderungen f schwa 
dieser Folie liefern zwei weitere Platinfolien mit großer Wandfläche. Aus den | des I 
Charakteristiken der Folien kann, da durch die Folien ein Strom bestimmter | haben 
Stärke fließen muß, wenn zwischen diese eine Spannung gelegt wird, die | [onen 
Änderung des Potentials der Absaugfolie gegenüber dem der Wand bei bekannter | werte 
Spannungsdifferenz bestimmt werden. Diese Potentialänderung ist in un- | auftre 
günstigsten Fällen beim Absaugen negativ geladener Teilchen bei kleinen Ent- f am gr 
ladungsstromstärken in der T-Form um 10% kleiner als die Spannungsdifferenz | In de 
zwischen den Folien. In der H-Form und bei der Untersuchung positiver | Werte 
Ionen entspricht sie praktisch der angelegten Spannung. die V 

Die Potentialschwankungen der Folien können die Messungen verfälschen, | meneı 
wenn zwischen den Schwingungsvorgängen auf der Absaugfolie und den | Verlaı 


Folien, die das Bezugspotential liefern, eine Phasenverschiebung auftritt. i 
Untersuchungen mit einem galvanisch gekoppelten Zweistrahl-Oszillographen | Energ 
ergaben, daß beide Vorgänge vollkommen in Phase sind. desser 


Da die Folien den Schwankungen des Plasmapotentials folgen, kann diese f werde 
Meßanordnung durch Anlegen eines Gegenfeldes zur Ermittlung der Energien | stets | 


der effundierenden Teilchen herangezogen werden. Entsprechende Messungen Da 
wurden durchgeführt. edi 
Ein regelbarer Hochspannungsgleichrichter liefert die Spannung fiir den ae 
Entladungsstromkreis, ein stabilisiertes Netzgerät die Beschleunigungsspannung | \erurs 
von 1800 Volt zwischen Effusionsöffnung und dem inneren Rohr der Rohr- a 
5) U. Weimer, Z. Naturforschg. 13a, 278 (1958). an Ne 


8) V.H. Attree, Energineering 21, 100 (1949) 
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linse und Anodenbatterien die Linsenspannung von 270 Volt und die außerdem 
noch benötigten Spannungen. Der Magnet des Massenspektrometers wird 
durch eine Akkumulatorenbatterie gespeist. 

Die Ionenströme am Auffänger werden durch eine von Kreuzer’) an- 
gegebene Brückenanordnung mit einer Elektrometerröhre gemessen. Die Nach- 
weisgrenze liegt bei 10-1 Amp. 


Die Meßgenauigkeit und Reproduzierbarkeit 


Die Meßgenauigkeit und Reproduzierbarkeit der Ergebnisse wird durch die 
Schwankungen der Entladungsparameter, durch die wechselnde Energie- 
breite der zu registrierenden Teilchen und durch die Konstanz der Trans- 
mission des Massenspektrometers festgelegt. 


Die Entladungsstromstärke, der Druck im Entladungsrohr und die Schwan- 
kungen des Plasmapotentials vor der Effusionsöffnung bestimmen, wie später 
gezeigt wird, wesentlich die Größe des Effusionsstromes. Stromstärke- 
schwankungen treten hauptsächlich in der 7’-Form auf, sie sind aber so unbe- 
deutend, daß sie die Meßergebnisse kaum beeinflussen dürften. Stärker 
machen sich Druckschwankungen bemerkbar. Obgleich sie nicht größer als 
0,1 Torr sind, können sie eine Änderung der Teilchenstréme um 5%, zur Folge 
haben. Diese Fehlerquelle kann nahezu ausgeschaltet werden, wenn nur un- 
mittelbar vor dem Gaseinlaß durch die Druckstabilisierungsanlage gemessen 
wird. 

Aus den MeBergebnissen kann geschlossen werden, daß die Potential- 
schwankungen des Plasmas vor der Austrittsöffnung am stärksten die Größe 
des Effusionsstromes beeinflussen. Die wechselnden Potentialverhältnisse 
haben zur Folge, daß zwischen verschiedenen Meßreihen bei den positiven 
[onen und Elektronen in ungünstigen Fällen Differenzen bis zu 30%, des Höchst- 
wertes und bei den negativen Ionen sogar solche bis zu 40% dieses Wertes 
auftreten können. Die Abweichungen sind bei allen Teilchen in der H-Form 
am größten, so daß die obigen Werte Höchstwerte in der H-Form angeben. 
In der T-Form sind die Abweichungen um 50%, niedriger. Wenn auch die 
Werte der Teilchenströme nicht vollkommen reproduzierbar sind, so sind doch 
die Verhältnisse der bei verschiedenen Entladungsstromstärken aufgenom- 
menen Teilchenströme stets nahezu die gleichen, so daß in ihrem qualitativen 
Verlauf keine wesentlichen Unterschiede auftreten. 

Eine mögliche Beeinträchtigung der Meßgenauigkeit durch die variable 
Energiebreite der Teilchen kann überprüft werden, indem der Gesamtstrom, 
dessen Wert auch bei mäßiger Fokussierung in seinem vollen Wert erfaßt 
werden kann, mit der Summe der lonenanteile verglichen wird. Es war 
stets Übereinstimmung vorhanden. 


Da nur bei einem Entladungsdruck von 2,5 Torr gemessen wurde, waren 
die durch das ausströmende Neutralgas bestimmten Druckverhältnisse im 
Massenspektrometer konstant und damit auch die durch gaskinetische Stöße 
verursachten Verluste. Der Absolutwert des Effusionsstromes wurde nicht 
gemessen, da er nicht interessierte. Er dürfte aber lediglich durch die Streuung 
an Neutralteilchen vermindert werden. denn durch die Potentialverteilung 


?) J. Kreuzer, Z. Physik 125, 707 (1949). ade. 
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in der Rohrlinse wird eine verlustfreie Bündelung des Teilchenstromes ge- 
währleistet. 

Bei erneutem Zusammenbau des Massenspektrometers konnte die Anlage 
so justiert werden, daß die Auffängerströme stets wieder in den durch die 
Schwankungen der Entladungsparameter bestimmten Grenzen lagen. 


Meßergebnisse 

Die Messungen wurden mit verschiedenen Effusionsöffnungen, die sich 
entweder in einer aufgeschmolzenen Platinfolie oder in der Rohrwand selbst 
befanden, durchgeführt. Die Durchmesser dieser Öffnungen schwankten 
zwischen 60 und 100 u. Da sich im Verlauf der Teilchenströme in Abhängig- 
keit von der Entladungsstromstärke bei Verwendung von Effusionsöffnungen 
in den verschiedenen Materialien keine wesentlichen Unterschiede zeigten, 
sind nur die Ströme, die mit einer Effusionsöffnung in einer aufgeschmolzenen 
Platinfolie ermittelt wurden, mitgeteilt. Der Durchmesser dieser Öffnung 
von 58 u muß als klein gegen die Dicke der Langmuirschicht angenommen 
werden. Nach dem von Pahl angegebenen Verfahren konnte hier aus den 
Meßwerten nicht auf die relativen Trägerkonzentrationen im Plasma ge- 
schlossen werden, da die Voraussetzungen der Schottkyschen Theorie, auf 
die sich die Pahlschen Überlegungen aufbauen, wegen der Anwesenheit von 
negativen Ionen im Plasma nicht erfüllt sind. Auch wird bei Sauerstoffglimm- 
entladungen der austretende Teilchenstrom stark durch die Potential- 
schwankungen des Plasmas beeinflußt, so daß die Diffusion zur Wand nicht 
allein die bestimmende Größe ist. Abb. 3 zeigt den Verlauf des Effusions- 
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20 40 60 80 700 J 
(mA) 
mL Abb. 3. Effusions- und Absaugstrom der Of-Ionen 

> 
stromes der O3-Ionen und den der extrahierten Träger in Abhängigkeit von 
der Entladungsstromstärke. Die Zahlen zwischen den Kurven geben die Diffe- 
renzen zwischen den Absaugspannungen gegenüber dem Wandpotential in 


Volt an. Auffallend ist die starke Abnahme des Ionenstromes beim Übergang 
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von der T- in die H-Form. Die MeBwerte und auch die folgenden wurden 
aufgenommen, nachdem die Entladung mit der entsprechenden Stromstärke 
eingebrannt worden war. Im wesentlichen war diese Maßnahme aber nur 
nach dem Sprung der Entladung von 14:0 

der T- in die H-Form erforderlich; 
denn hier trat während des Einbrennens 


4 
sowohl bei den positiven Ionen als 
auch bei anderen Trägern eine merkliche _,, | 
Abnahme des Teilchenstromes auf, so 
daß erst nach einigen Minuten der 0 
Endwert erreicht wurde. Die massen- ] 
spektrometrische Untersuchung der 
positiven Teilchenströme ergab, daß in 67 
der T-Form nur O3-Ionen effundieren 
und extrahiert werden können. Weitere 61 


Ionenarten konnten nicht nachgewiesen 
werden, obgleich die Empfindlichkeit 44 
und das Auflösungsvermögen der An- 
ordnung es erlaubte, Teilchenströme 24 
anderer Ionen .in der Größenordnung 
von 10-13 Amp. zu registrieren. Falls -— 
noch andere lIonenarten auftreten, 80 100 20 J 
muß das Verhältnis ihrer Ströme zum 
Molekülionenstrom kleiner als 1:1000 
sein. 

In der H-Form bestimmen ebenfalls die OF-Ionen wesentlich den Ionen- 
strom. Es konnte daneben noch eine schwache O{-Komponente gemessen 
werden, deren Größe in Abhängigkeit von der Entladungsstromstärke und der 
Ziehspannung in Abb. 4 angegeben ist. 

Der Elektronenstrom (Abb. 5) zeigt praktisch den gleichen Verlauf wie 
der positive Ionenstrom; nur ist sein Wert häufig größer. Maximal betrug die 


Abb. 4. Effusions- und Absaugstrom der 
O#+-Ionen 


20 40 60 80 100 
0 Abb. 5, Effusions- und Absaugstrom der Elektronen 
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Differenz 50%, des totalen Ionenstromes. Diese Abweichung beens PR 
hervorgerufen werden, daß durch eine geringe Deformation der Feldverteilung 
in der Langmuirschicht, die auf den Ionenstrahl zerstreuend wirkt, die Ver- 
luste der Ionen im Effusionskanal größer sind als die der Elektronen. 


Abb. 6. Effusions- und Absaugstrom der O--Tonen 


eh 


80 100 20 J 


Abb. 7. Effusions- und Absaugstrom der O, -lonen 


Negative Ionen konnten sowohl in der T-Form als auch in der H-Form 
nachgewiesen werden. Es wurden O-- und O,-Ionen in beiden Entladungs- 
formen registriert, wobei die O,-Ionen in der H-Form nur bei Extraktion mit 
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einer Ziehspannung von 2 Volt gefunden wurden. Im Verlauf der hier ange- 
gebenen Teilchenströme (Abb. 6 und 7) treten Abweichungen von dem der 
positiven Ionen und Elektronen auf. Zwar ist der Sprung beim Übergang 
von der 7'- in die H-Form genauso vorhanden; aber bei niedrigen Entladungs- 
stromstärken sind Unterschiede in der 7'-Form bemerkbar. Bei beiden Ionen- 
arten ist bei einer Stromstärke von 20 mA ein Maximum erkennbar. Neben 
den hier gezeigten Kurven mit den stark ausgeprägten Maxima konnten aber 
auch solche aufgenommen werden, bei denen das Maximum nur schwer bzw. 
gar nicht erkennbar war, und zwar dann, wenn das Sprunggebiet unter sehr 
sauberen Verhältnissen bei hohen Stromstärken auftrat. 


Die Energien der austretenden Teilchen i 


Durch Gegenfeldmessungen konnte die Energie der effundierenden Ionen 
und Elektronen ermittelt werden. Die N werden normiert in 


ist die groBe Differenz mA 
energien in der H- und ser Jah 
T-Form und die Paralleli- 
energie und der Ampli- 
tude der Potentialschwan- 15 2 u 

kungen (Abb. 12). iv 


Die Potentialschwan- Abb.8. Energie der effundierenden positiven Ionen 
kungen, die mit galva- 

nisch gekoppelten Oszillographen gemessen wurden, treten in der 7'-Form 
teilweise als Schwingung auf, denen statistische Potentialschwankungen über- 


lagert sind. In der H-Form sind die Verhältnisse übersichtlicher; denn es bilden 
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Lange Abb. 9. Energie der effundierenden Elektronen 


sich in dem hier untersuchten Bereich reine Schwingungsvorgänge aus. In 
der T-Form ließen sich die Amplituden wegen der unübersichtlichen Potential- 
schwankung schwer bestimmen. Es sind hier die Werte der erkennbaren 
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Schwingungsvorgänge, nicht aber die Spitzenpotentiale der überlagerten 
statistischen Schwankungen angegeben. Der in Abb. 12 durch die ausgezogene 


Kurve angegebene Amplitudenverlauf wurde unter Verhältnissen aufgenom- 
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Abb. 11. Energie der effundierenden O--Ionen 
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: men, bei denen das Sprunggebiet bei Entladungsstromstärken liegt, die höher 

7a als 80 mA sind, also unter sauberen Entladungsbedingungen. Es war jedoch 
nicht möglich, derartig saubere Verhältnisse für alle Versuchsreihen herzu- 
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stellen. Treten geringe Verunreinigungen auf, dann verschiebt sich das 
Sprunggebiet zu kleineren Stromstärkewerten und bei einer Entladungs- 
stromstärke von 20 bis 30 mA bildet sich im Potentialverlauf ein Maximum 
aus. Der Potentialanstieg kurz vor dem Übergang von der T- in die H-Form 
ist dann ebenfalls stärker ausgeprägt. Dieser Verlauf ist in Abb. 12 gestrichelt 
angedeutet. 


Die mit der Gegenfeldmethode gemessenen Ionenenergien sind innerhalb 
der MeBgeriauigkeit sowohl im Wert als auch im Verlauf in guter Uberein- 
stimmung mit der Größe der Amplitude der Potentialschwingungen vor der 
Effusionsöffnung. Es konnte in der H-Form, wo die Verhältnisse am über- 
sichtlichsten sind, gezeigt werden, daß strenge Proportionalität zwischen der 
Amplitude der Potentialschwingungen und der Teilchenenergie besteht. 
Es verdoppelt sich kurz vor dem Sprung von der H- in die T-Form die Ampli- 
tude der Schwingungen und ebenfalls auch die Energie der austretenden 
Ionen. 


Da die Absaugfolie und auch die Folien, die das Bezugspotential liefern, 
den Potentialschwankungen wie unbelastete Folien folgen und außerdem 
zwischen den Schwingungsvorgängen auf den entsprechenden Folien keine 
Phasenverschiebung besteht, so folgt aus den Meßwerten, daß entweder die 
Potentialschwingungen vor der Effusionsöffnung senkrecht zur Wand auf- 
treten, oder aber daß die Folien mit einer Phasenverschiebung von 180° den 
Potentialschwingungen des Plasmas folgen. Da die letzte Annahme sehr un- 
wahrscheinlich ist, so darf man vermuten, daß die Potentialschwingungen in 
radialer Richtung erfolgen. Völlig sicher kann man diese Schwingungs- 
komponente aber nur mit einem Elektronenstrahl nachweisen, so wie es 
Gabor, Ash und Dracott®) am Quecksilber-Niederdruckbogen durchführten. 
Sie konnten mit einem Elektronenstrahl, den sie parallel zur Wand aber 
senkrecht zur Entladungsachse durch den interessierenden Bereich schickten, 
zeigen, daß sich in solchen Entladungen Potentialschwingungen mit einer 
Komponente senkrecht zur Wand ausbilden. 


Bei der Untersuchung der Effusionsströme verändert sich die Amplitude 
der Potentialschwankungen mit der Entladungsstromstärke und entsprechend 
auch die Energie. der effundierenden Teilchen; bei positiven Ionen aller- 
dings nicht in so starkem Maße wie die Amplitude der Potentialschwankungen. 
Die Abweichungen liegen innerhalb der Meßgenauigkeit; denn der Absolut- 
wert der Ionenenergie kann mit einem Fehler von 30%, behaftet sein, während 
bei den Energiedifferenzen der verschiedenen Ionenarten Abweichungen bis 
zu 100%, auftreten können. Die in der H-Form nachgewiesenen Of-Ionen 
haben stets eine höhere Energie als die O$-Ionen. Die Differenz beträgt im 
Durchschnitt 0,5eV. Sie kann durch die Unterschiede der Beweglichkeits- 
werte bedingt sein. Burch und Geballe®) konnten für O$-Ionen, deren Be- 
wegung durch quantenmechanische Effekte beeinflußt wird, einen Beweglich- 
keitswert von 1,95 cm?/Vsec messen. Für Of-Ionen berechneten sie nach 
der Langevinschen Theorie, ohne Ladungsaustausch und ae 
effekte berücksichtigen zu müssen, einen Wert von 2,34 cm?/Vsec. is ; 


8) D. Gabor, E. A. Ash u. Dr. Dracott, Nature 176, 916 (1955). 
®) D.S. Burch; R. Geballe, Phys. Rev. 106, 183 (1957). 
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Die Energie der austretenden Elektronen befindet sich in guter — ind 1 
stimmung mit der Amplitude der Potentialschwankungen. ot i 
Burch und Geballe®) und andere Autoren zeigten, daß praktisch keine = 
Differenz zwischen den Beweglichkeitswerten der entsprechenden positiven 
und negativen Ionen besteht. Sie gaben für O--Ionen einen Meßwert von De 
3,4 em?/Vsec und für die O2-Ionen einen solchen von 1,95 em?/Vsec an. Diese wert, 
Differenz macht sich auch bei den Energien der entsprechenden Ionen be. | @"- I 
merkbar, und zwar ist die Energie der O2-Ionen stets kleiner als die der O-. dieren 
Ionen. Der Unterschied hat im Durchschnitt eine Größe von | eV. Die bei Di 
verschiedenen Entladungsstromstärken gemessenen Energien der negativen Größe 
Ionen stehen qualitativ in guter Übereinstimmung mit der Amplitude der | Entla 
Potentialschwankungen. tentia 
Der Effusionsstrom der positiven und negativen Träger wird durch die gestal 
Schwingungsamplitude der Potentialschwankungen wesentlich beeinflußt; Elekt 
denn neben dem qualitativ gleichen Verlauf beider Größen darf man die | stärke 
Tatsache als Beweis hierfür werten, daß, wenn sich kurz vor dem Sprung | Entla 
von der H- in die T-Form die Amplitude der Potentialschwingungen ver- kleine 
doppelt, der Effusionsstrom der positiven Ionen und Elektronen sich eben- 
falls erhöht, auch wenn die Effusionsöffnung sich in einer Glaswand be- 
findet. Es ist daher der Effusionsstrom nicht unmittelbar ein Maß für die A 
Trägerkonzentration in der Nähe der Rohrwand. 
vad ae tratic 
Der Absaugstrom 
Mit einer Ziehspannung bis zu 20 V gegenüber dem Wandpotential wurden auch 
geladene Träger extrahiert. Der Verlauf der Absaugströme in der T- = H- ward 
Form zeigt wesentliche Unterschiede. (Abb. 13 und 14.) Sis aah. wee M 
weite 
/ 36 bis zı 
3. Absaugstrom der posi- Abb. 14. Absaugstrom der Stroi 
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Völlig gleichartig ist das Verhalten der negativ geladenen Teilchen, darum 
sind hier nur die Ströme der Elektronen mitgeteilt. Der Absaugstrom der 
O-Ionen zeigt weder in der H- noch in der 7T-Form einen Sättigungswert, 
während bei dem der Of-Ionen Sättigungserscheinungen angedeutet sind. 

Der Strom der negativen Teilchen erreicht in der T-Form einen Sättigungs- 
wert, in der H-Form steigt er in einem Spannungsbereich von 12—20 V steil 
an. Die Absaugströme der Teilchen, die in beiden Entladungsformen effun- 
dieren, haben in der H-Form eine stärkere Zunahme als in der T-Form. 


Die Absaugkurven können hier nicht näher diskutiert werden, da die 
Größen, die die Wandströme bestimmen, in ihren Werten in den jeweiligen 
Entladungsformen nicht bekannt sind. Außerdem läßt sich wegen der Po- 
tentialschwankungen, die die Verhältnisse vor der Rohrwand unübersichtlich 
gestalten, schwer angeben, durch welche Gesetzmäßigkeiten die Ströme der 
Elektronen und Ionen zur Wand geregelt werden. Man darf aber wegen des 
stärkeren Anstiegs der Absaugströme in der H-Form vermuten, daß in dieser 
Entladungsform die Trägerkonzentration vor der Effusionsöffnung nicht 
kleiner ist als die in der 7'-Form. 


Schlußbemerkungen 


Aus den Werten der Effusionsströme kann nicht auf die Trägerkonzen- 
trationen im Plasma geschlossen werden, da keine Aussagen iiber die Zu- 
sammenhänge zwischen den Effusionsströmen und den relativen Träger- 
konzentrationen vorliegen, die auf Grund theoretischer Überlegungen, die 
auch die Schwingungsvorgänge in der Entladung berücksichtigen, gemacht 
wurden. 


Man darf annehmen, daß neben den hier gefundenen Ionenarten noch 
weitere im Entladungsplasma auftreten. So sind neben den O%-Ionen noch 
O+-Ionen zu erwarten, denn spektroskopische Untersuchungen zeigten neben 
den O,-Banden auch O-Linien!). Welche Vorgänge die Effusion der O+-Ionen 
verhindern, konnte nicht geklärt werden. Es ist unwahrscheinlich, daß die 
Schwingungserscheinung vor der Rohrwand als Ionenfilter wirkt; denn die 
wechselnde Frequenz der Potentialschwingungen, die sich von einigen kHz 
bis zu 40 kHz erstreckt, läßt eine gleichbleibende Filterung als unmöglich er- 
scheinen. Außerdem widerlegt die Effusion der O--Ionen diese Annahme. 


Der Verlauf der Effusionsströme und der Amplitude der Plasmaschwin- 
gungen läßt erwarten, daß die Trägerkonzentration an O3-, O--Ionen und 
Elektronen in der T-Form proportional mit der Entladungsstromstärke an- 
steigt, wobei bei dem Umschlag der Entladungsform keine Abnahme der 
Trägerkonzentration auftritt. In der H-Form dürfte die Größe des Effusions- 
stromes ein Maß für die relativen Trägerkonzentrationen an der Rohrwand 
sein, da die Schwingungsamplitude praktisch konstant bleibt. Wesentliche 
Verschiedenheiten treten in der H-Form im Verhalten der O$- und der Of-Ionen 
auf. Während der O3-Strom praktisch unabhängig von der Entladungsstrom- 
stärke ist, erhöht sich die Zahl der effundierenden Of-Ionen, außerdem läßt 
der Verlauf des Absaugstromes auf einen Sättigungswert schließen. 

Die Konzentration der O -Ionen nimmt in der H-Form mit steigender 
Stromstärke ab, während die der O>-Ionen konstant bleibt. Bei dem Uber- 
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gang von der T- in die H-Form kann eine sprunghafte Geier der 
O2-Konzentration als sicher angenommen werden. 


Herrn Professor Dr. Seeliger möchte ich für die Anregung zu dieser 
Arbeit und sein Interesse an ihrem Fortgang herzlich danken. Herrn Dr. 
Wulfhekel bin ich für sein stetes Interesse sehr zu Dank verpflichtet. 


Greifswald, Institut für Gasentladungsphysik der Deutschen Akademie 


der Wissenschaften zu Berlin. 


A ye Bei der Redaktion eingegangen am 17. August 1959. gi ste 
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HORT-THOMA 


7. Auflage des Lehrbuches 
„Die Differentialgleichungen des Ingenieurs“ von W. HORT 


Neubearbeitet von Oberstudiendirektor Dr. phil. ALFRED THOMA, Fulda 
XI, 582 Seiten mit 343 Abb. im Text. 1956. gr. 8°. DM 36.—, Leinen. DM 38.— 


„Zentralblatt für Mathematik und ihre Grenzgebiete‘: ,,Der Verfasser hat mit 
viel Geschick und Umsicht die Fülle des den verschiedenen Gebieten der Technik und 
Physik entnommenen Stoffes zu einem Ganzen geformt, ohne dabei das Ziel aus den Augen 
zu verlieren, den Nutzen der Mathematik im Schaffen des Ingenieurs deutlich erkennbar 
zu machen. Da das Buch in erster Linie für Ingenieure und Physiker geschrieben ist, liegt 
bei der Darstellung der mathematischen Methoden die Betonung auf dem Handwerklichen. 
Die vorliegende Auflage des bewährten Werkes wird sich sicherlich in Ingenieurkreisen 
viele Freunde erwerben.‘ W. Quade 


Operatorenrechnung 


Laplacesche Transformation nebst Anwendungen 


in Physik und Technik 
Von Prof. Dr. phil., Dr.-Ing. e. h. K. W. WAGNERT 
2., verbesserte Auflage 


XVIII, 471 Seiten mit 126 Abb. im Text und einem Porträt von Laplace. 
1950. gr. 8°. DM 40.30, Leinen DM 42.80 


VDI-Zeitschrift: Der Verfasser führt in die ganze Theorie ein, indem er sehr 
wertvolle praktische Beispiele aus verschiedenen Disziplinen der Naturwissenschaft und 
Technik durchrechnet. Man lernt, die Verfahren zu handhaben, aus dem Gefühl heraus, 
daß sie die einzige Möglichkeit bieten, die Probleme anzugreifen. Das Studium des vor- 
liegenden Werkes kann Naturwissenschaftlern, Technikern und Mathematikern wärmstens 
empfohlen werden. U, Wegner 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 
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Theoretische Chemie 


a Eine Einführung vom Standpunkt einer gestalthaften Atomlehre 
Von Prof. Dr. K. LOTHAR WOLF, Kirchheimbolanden 

4., umgearbeitete Auflage 
XI, 813 Seiten mit 347 Abb. im Text. 1959. gr. 8°. DM 39.50. Leinen DM 41.50 


Chemiker-Zeitung: Das Werk, das sich bereits jetzt einen festen Platz in der 
physikalisch-chemischen Literatur erobert hat, schreitet weiter auf seinem eigenwilligen 
Weg, der bewußt ein anderer ist als der, den die ‚theoretische Chemie thermo-dyna- 
mischer Prägung‘ geht. Das Werk bietet handfeste physikalische Chemie, und der 
Leser — auch Anfängern kann es empfohlen werden —, der es durcharbeitet, lernt nicht 
nur viel theoretisch-chemischen Stoff, sondern er darf auch die Freude an einer Dar- 
stellungsweise genießen, die ihm den Standort der Einzelerscheinungen besonders klar 
aufzeigt. Dem Werk ist eine weite Verbreitung zu wünschen. Sch. 


Praktische Einführung 
In die physikalische Chemie 
Von Prof. Dr. K. LOTHAR WOLF, Kirchheimbolanden 
und Dr. hatil. HANS GEORG TRIESCHMANN, Ludwigshafen/Rhein 


3., neubearbeitete Auflage 
VIII, 255 Seiten mit 66 Abb. i. Text. 1954. 8°. Halbleinen DM 9.70 


Chemiker-Zeitung: Dasin dritter Auflage vorliegende kleine Werk ist nicht nur zum 
Gebrauch beim physikalisch-chemischen Praktikum, sondern auch — und dies ist sein 
besonderes Anliegen — als Hilfsmittel zu lebendiger Veranschaulichung des Lehrstoffes 
bei der Durcharbeitung größerer Lehrbücher der physikalischen Chemie und neben det 
Vorlesung sowie als selbständige kurze Einführung in die physikalische Chemie gut 


geeignet Sch. 


Kristallchemie und Kristallphysik 
metallischer Werkstoffe 


Von Dr. FRANZ HALLA, Briissel 
Em. Prof. an der Technischen Hochschule Wien 
Wissenschaftlicher Mitarbeiter der Association pour les Etudes Texturales, Brüssel 


3., völlig umgearbeitete Auflage 
XVI, 737 Seiten mit 365 Abbildungen im Text und auf einer Tafel. 1957. gr. 8°. 
DM 50.—, Leinen DM 52.— 
Angewandte Chemie: Das Werk stellt eine bewundernswerte Leistung des Verfassers 
dar, der aus zahlreichen Einzelveröffentlichungen ein außerordentlich großes Tatsachen- 
material zusammengetragen hat. Das Buch bildet ein wertvolles Bindeglied zwischen 
voneinander weit entfernten Arbeitsgebieten der Physik, Chemie, Metallkunde und der 


Technik. Es kommt daher für einen großen Interessentenkreis in Frage und wird als 
Lehr- und Nachschlagewerk mit Nutzen verwendet werden. G. Wassermann 


JOHANN AMBROSIUS BARTH / VERLAG / LEIPZIG 
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